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THEORETICAL  ANALYSIS  OF  MICROWAVE  AND  MILLIMETER  WAVE  INTEGRATED 


CIRCUITS  BASED  ON  MAGNETIC  FILMS 


Under  the  sponsorship  of  the  ONR  Contract  Contract  N00014-89-J-1019  we  have 
published  24  refereed  journal  and  conference  papers. 

A  full  modal  analysis  is  used  to  study  the  dispersion  characteristics  of  microstrip 
lines  periodically  loaded  with  crossing  strips  in  a  stratified  uniajciadly  anisotropic  medium. 
Dyadic  Green’s  functions  in  the  spectral  domain  for  the  multilayered  medium  in  conjunc¬ 
tion  with  the  vector  Fourier  transform  (VFT)  are  used  to  formulate  a  coupled  set  of  vector 
integral  equations  for  the  current  distribution  on  the  signal  line  and  the  crossing  strips. 
Galerkin’s  procedure  is  applied  to  derive  the  eigenvalue  equation  for  the  propagation  con¬ 
stant.  The  effect  of  anisotropy  for  both  open  and  shielded  stuctures  on  the  stopband 
properties  is  investigated. 

The  input  imped2Lnce  of  a  microstrip  antenna  consisting  of  two  circular  microstrip 
disks  in  a  stacked  configuration  driven  by  a  coaxial  probe  is  investigated.  A  rigorous  anal¬ 
ysis  is  performed  using  a  dyadic  Green’s  function  formulation  where  the  mixed  boundary 
value  problem  is  reduced  to  a  set  of  coupled  vector  integral  equations  using  the  vector 
Hankel  transform.  Galerkin’s  method  is  employed  in  the  spectral  domain  where  two  sets 
of  disk  current  expansions  are  used.  One  set  is  based  on  the  complete  set  of  orthogonal 
modes  of  the  magnetic  cavity,  and  the  other  employs  Chebyshev  polynomials  with  the 
proper  edge  condition  for  the  disk  currents.  An  additional  term  is  added  to  the  disk  cur¬ 
rent  expansion  to  properly  model  the  current  in  the  vicinity  of  the  probe/disk  junction. 
The  input  impedance  of  the  stacked  microstrip  antenna  including  the  probe  self-impedance 
is  calculated  as  a  function  of  the  layered  parameters  and  the  ratio  of  the  two  disk  radii. 
Disk  current  distributions  and  radiation  patterns  are  also  presented.  The  calculate  results 
are  compared  with  experimentEd  data  and  shown  to  be  in  good  agreement. 
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The  coupled-wave  theory  is  generalized  to  analyze  the  diffraction  of  waves  by  chiral 
gratings  for  arbitrary  angles  of  incidence  and  polcirizations.  Numerical  results  for  the 
Stokes  peirameters  of  diffracted  Floquet  modes  versus  the  thickness  of  chiral  gratings  with 
various  chiralities  are  calculated.  Both  horizontal  and  vertical  incidences  are  considered  for 
illustration.  The  diffracted  waves  from  chiral  gratings  are  in  general  elliptically  polarized; 
and  in  some  particular  instances,  it  is  possible  for  chiral  gratings  to  convert  a  linesu’ly 
poleirized  incident  field  into  two  nearly  circularly  polcirized  Floquet  modes  propagating  in 
different  directions. 

A  general  spectral  domain  formiilation  to  the  problem  of  radiation  of  arbitrary 
distribution  of  sources  embedded  in  a  horizontally  stratified  arbitrary  magnetized  linear 
plasma  is  presented.  The  fields  are  obtained  in  terms  of  electric  and  magnetic  type  dyadic 
Green’s  functions.  The  formulation  is  considerably  simplified  by  using  the  kDB  system 
of  coordinates  in  conjunction  with  the  Fourier  transform.  The  distributional  singul^u: 
behavior  of  the  various  dyadic  Green’s  functions  in  the  source  region  is  investigated  and 
taken  into  account  by  extracting  the  delta  function  singidarities.  Finally,  the  fields  in  any 
arbitrary  layer  are  obtained  in  terms  of  appropriately  defined  globed  upward  and  downward 
reflection  and  transmission  matrices. 

We  have  investigated  a  method  for  the  calculation  of  the  current  distribution,  resis¬ 
tance,  and  inductance  matrices  for  a  system  of  coupled  superconducting  transmission  lines 
having  finite  rectangular  cross  section.  These  cedculation  allow  accurate  characterization 
of  both  high-Te  and  low-T*  superconducting  strip  transmission  lines.  For  a  single  stripline 
geometry  with  finite  ground  planes,  the  current  distribution,  resistance,  inductance,  and 
kinetic  inductance  are  calciilated  as  a  function  of  the  penetration  depth  for  various  film 
thickness.  These  cedculations  are  then  used  to  determine  the  penetration  depth  for  Nb,  NbN, 
and  Y BaiCuzOj-t  superconducting  thin  films  from  the  measured  temperature  dependence 
of  the  resonant  frequency  of  a  stripline  resonator.  The  calculations  are  also  used  to  convert 
measured  temperature  dependence  of  the  quality  factor  to  the  intrinsic  surface  resistance 
as  a  function  of  temperature  for  a  Nb  stripline  resonator. 
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Proximity-print  x-ray  lithography  is  commonly  performed  with  gold  or  tungsten 
structures  of  sizes  down  to  30  nm  wide  and  50-800  nm  tall  which  are  patterned  onto  the 
surface  of  a  thin,  x-ray  transparent  membrane.  X-rays  in  the  wavelength  range  of  0.5-5  nm 
are  used  for  replication  with  mask-substrate  gaps  ranging  £rom  zero  (contact  print)  up  to 
20  nm  or  more.  The  resolution  o  f  this  method  (minimum  achievable  linewidth)  is  limited 
predominantly  by  the  diffraction  of  the  x-rays  around  these  structures  and  the  spreading 
of  the  diffracted  waves  into  the  0—20  gap.  Work  to  date  hjis  assumed  that  scalar 
diffraction  theory  is  applicable — as  calculated,  for  example,  by  the  Rayleigh- Sommerfeld 
formulation — and  that  Kirchhoff  boundary  conditions  can  be  applied.  Kirchhoff  boundary 
conditions  assume  that  the  fields  are  constant  in  the  region  between  the  absorbers,  and 
also  (a  different)  constant  in  the  region  just  under  the  absorbers,  and  that  there  are  no 
fringing  fields.  In  this  report  we  explore  the  validity  of  this  assumption  for  the  case 
of  30  nm-wide  by  30-100  nm-tall  gold  absorbers  with  4.5  nm  {Ck)  x-rays.  Because  of 
computational  time  limitations,  the  shorter  wavelength  and  larger  absorber  cases  are  not 
currently  possible.)  Because  the  absorber  is  only  7  wavelengths  wide  and  7-20  wavelengths 
high,  strong  diffractive  effects  are  expected.  The  finite-difference  time-domain  (FD-TD) 
technique  was  used  on  a  Cray-2  supercomputer  to  predict  the  fields  diffracted  by  the  gold 
absorbers.  In  applying  the  FD-TD  technique.  Maxwell’s  equations  are  discretized  in  space 
and  time  on  a  uniform  rectangular  grid.  A  second-order  absorbing  boundary  condition  is 
applied  at  the  outer  boundary  of  the  computational  domain  in  order  to  simulate  unbounded 
space.  The  results  indicate  that  strong  fringing  fields  exist  in  the  shadow  region  of  the 
absorber,  and  hence  Kirchhoff  boundary  conditions  are  not  accurate  in  this  regime. 

Because  the  effects  of  diffraction  during  proximity-print  x-ray  lithography  are  of 
critical  importance,  a  number  of  previous  researchers  have  attempted  to  calculate  the 
diffraction  patterns  and  minimum  achievable  feature  sizes  as  a  function  of  wavelength  and 
gap.  Work  to  date  has  assumed  that  scalar  diffraction  theory  is  applicable— as  calculated,  for 
example,  by  the  Rayleigh-Sommerfeldformulation-and  that  Kirchhoff  boundary  conditions 
can  be  applied.  Kirchhoff  boundary  conditions  assume  that  the  fields  (amplitude  and 
phase)  are  constant  in  the  open  regions  between  absorbers,  and  a  different  constant  in 
regions  just  under  the  absorbers  (i.e.,  that  there  are  no  fringing  fields).  An  x-ray  absorber 
is,  however,  best  described  as  a  lossy  dielectric  that  is  tens  or  hundreds  of  wavelengths 
tall,  and  hence  Kirchhoff  boundary  conditions  are  unsuitable.  In  this  report  we  use  two 
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numerical  techniques  to  c^dculate  (on  a  Cray  2  supercomputer)  accurate  diffracted  fields 
from  gold  absorbers  for  two  cases:  a  30  nm-wide  line  at  A  =  4.5  nm,  and  a  100  nm-wide 
line  at  A  =  1.3  nm.  We  show  that  the  use  of  Kirchhoff  boundary  conditions  introduces 
unphysically  high  spatial  frequencies  into  the  diffracted  fields.  The  suppression  of  these 
frequencies- which  occurs  naturally  without  the  need  to  introduce  an  extended  source  or 
broad  spectrum-improves  exposure  latitude  for  mask  features  near  0.1  pm  and  below. 

In  order  to  understand  the  physical  meaning  of  rational  reflection  coefficients  in 
one-dimensional  inverse  scattering  theory  for  optical  waveguide  design,  we  have  studied 
the  relation  between  the  poles  of  the  transverse  reflection  coefficient  and  the  modes  in 
inhomogeneous  dielectrics.  By  using  a  stratified  medium  model  it  is  shown  that  these 
poles  of  the  reflection  coefficient  have  a  one-to-one  correspondence  to  the  discrete  modes, 
which  are  the  guided  and  leaky  modes.  The  radiation  modes  have  continuous  real  values  of 
transverse  wave  numbers  and  are  not  represented  by  the  poles  of  the  reflection  coefficient. 
Based  on  these  results,  applications  of  the  Gel’fand-Levit2in-Marchenko  theory  to  optical 
waveguide  synthesis  with  the  rational  function  representation  of  the  transverse  reflection 
coefficient  are  discussed. 

We  developed  an  inversion  algorithm  based  on  a  recently  developed  inversion  method 
referred  to  as  the  Renormalized  Source-Type  IntegreJ  Equation  approach.  The  objective 
of  this  method  is  to  overcome  some  of  the  limitations  and  difficulties  of  the  iterative  Born 
technique.  It  recasts  the  inversion,  which  is  nonlinear  in  nature,  in  terms  of  the  solution 
of  a  set  of  linear  equations;  however,  the  final  inversion  equation  is  still  nonlineM.  The 
derived  inversion  equation  is  an  exact  equation  which  sums  up  the  iterative  Neuman  (or 
Born)  series  in  a  closed  form  and;  thus,  is  a  valid  representation  even  in  the  case  when 
the  Born  series  diverges;  hence,  the  name  Renormalized  Source-Type  Integral  Equation 
Approach. 

The  scattering  and  receiving  characteristics  of  a  probe-fed  stacked  circular  mi¬ 
crostrip  antenna,  both  as  an  isolated  element  and  in  an  infinite  array,  are  investigated. 
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The  receiving  case,  where  the  antenna  is  loaded  with  impedance  is  solved  by 
superposition,  decomposing  the  problem  into  the  scattering  case  with  =  0  and  the  trans¬ 
mitting  case.  In  the  scattering  case,  the  coaxial  probe  is  short-circuited  to  the  ground  plane 
and  the  induced  probe  current  h  due  to  an  incident  plane  wave  excitation  is  determined. 
In  the  transmitting  case,  a  voltage  V  is  applied  to  the  base  of  the  probe  and  the  input 
impedance  is  solved  for,  giving  a  relationship  between  the  applied  voltage  V  and  the 
transmitting  probe  current  7j.  With  the  knowledge  of  h  and  Za,,  for  a  given  load  impedance 
Zl,  the  total  probe  current,  I  =  Ii+h,  and  the  received  power  are  determined. 

The  scattering  and  transmitting  problems  are  solved  rigorously  using  a  dyadic 
Green’s  function  formulation  where  the  mixed  boundary  value  problem  is  reduced  to  a  set 
of  coupled  vector  integral  equations  for  the  unknown  disk  and  probe  currents.  Galerkin’s 
method  is  employed  in  the  spectral  domain  where  the  disk  current  distributions  are  ex¬ 
panded  in  terms  of  the  complete  set  of  transverse  magnetic  (TM)  and  transverse  electric 
(TE)  modes  of  a  cylindrical  resonant  cavity  with  magnetic  side  wzdls.  An  addition2J  term 
is  added  to  the  disk  current  expansion  to  properly  model  the  sing\ilar  behavior  of  the 
current  in  the  vicinity  of  the  probe,  to  ensure  continuity  of  the  current  at  the  probe/disk 
junction,  and  to  speed  up  the  convergence  of  the  solution. 

The  radar  cross  section  (RCS)  of  a  single  stacked  microstrip  antenna  is  calculated 
for  both  the  open  and  short-circuited  cases.  For  an  infinite  array  of  phased  elements,  the 
reflection  coefficient  seen  at  the  input  of  the  antenna  and  the  received  power  are  calculated. 

The  complex  resonant  frequencies  of  the  open  structure  v:>f  a  microstrip  antenna 
consisting  of  two  circular  microstrip  disks  in  a  three  layer  stacked  configuration  have  been 
rigorously  calculated  as  a  function  of  the  layered  parameters  and  the  ratio  of  the  radii 
of  the  two  disks.  Using  a  dyadic  Green’s  function  formulation  for  horizontally  stratified 
media  and  the  vector  Hankel  transform,  the  mixed  boundary  value  problem  is  reduced  to 
a  set  of  coupled  vector  integral  equations.  Employing  Galerkin’s  method  in  the  spectral 
domain,  the  complex  resonant  frequencies  are  calculated  and  convergence  of  the  results  is 
demonstrated.  It  is  shown  that  for  each  mode,  the  stacked  circular  microstrip  structure  has 
dual  resonant  frequencies  which  are  cissociated  with  the  two  coupled  constitutive  resonators 
of  the  structure  and  which  are  a  function  of  the  mutual  coupling  between  them.  This  mutal 
coupling  depends  on  the  geometrical  configuration  of  the  stacked  structure,  the  layered 
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parameters,  and  the  disk  radii.  The  maximum  coupling  effect  occurs  where  the  real  parts 
of  the  resonant  frequencies  of  the  constitutive  resonators  are  approximately  equal,  where 
the  behavior  of  the  resonances  in  this  region  is  a  function  of  the  coupling.  The  dual 
frequency  behavior  of  the  stacked  microstrip  structure,  easily  controlled  by  varying  the 
parameters  of  layer  2  and  disk  radii  ratio,  given  fixed  parameters  for  layer  1  and  layer  3, 
may  be  used  to  broaden  the  bandwidth  or  provide  for  dual  frequency  use  of  the  antenna. 

We  rigorously  analyze  the  radiation  problem  of  a  circular  patch  which  is  center  fed 
by  a  coajdal-line  driven  probe  over  a  ground  plame  and  situated  in  an  arbitrary  layered 
medium.  The  current  distribution  on  both  the  patch  and  the  probe  is  rigorously  formu¬ 
lated  using  a  planar  stratified  medium  approach.  A  set  of  three  coupled  integral  equation 
is  derived  which  governs  the  axial  current  distribution  on  the  probe,  the  radial  current 
distribution  on  the  patch  and  the  sizimuthal  magnetic  current  sheet  across  the  aperture  of 
the  driving  coaxial  line.  This  set  of  equations  is  then  solved  using  the  method  of  moments. 
The  resulting  matrix  equation  is  obtained  in  terms  of  Sommerfeld-type  integrals  that  take 
into  account  the  effect  of  the  layered  medium.  These  integrals  are  efficiently  computed  by 
a  simple  deformation  in  the  complex  wavenumber  domain.  The  probe  current  distribu¬ 
tion,  input  impedance  and  radiation  pattern  are  presented  and  compaxed  to  the  case  of  a 
uniform  probe  current  distribution. 

Microstrip  antennas  of  stacked  configurations  have  received  attention  in  recent  years 
for  both  wideband  and  dual  frequency  use,  overcoming  the  narrow  bandwidth  of  conven¬ 
tional  single  layer  microstrip  antennas.  Although  much  experimental  work  has  been  per¬ 
formed,  theoretical  analyses  of  stacked  microstrip  patches  is  limited.  Resonant  frequencies 
of  the  stacked  microstrip  antennas  have  been  rigorously  calculated.  Numerical  methods 
have  been  used  to  calculate  the  current  and  radiation  fields  of  a  stacked  microstrip  antenna. 
The  method  of  moments  has  been  applied  to  analyze  to  the  stacked  microstrip  structure 
when  excited  by  an  incident  plane  wave.  A  spectral  domain  iterative  analysis  for  a  stacked 
microstrip  antenna  where  the  antenna  is  described  by  a  rectangular  sampling  grid  has  been 
used  to  calculate  radiation  patterns.  This  analysis  does  not  allow  for  accurate  modeling 
of  the  probe  feed. 


7 


In  our  approach,  the  input  impedance  and  radiation  fields  of  a  coaxial  probe-fed 
microstrip  antenna  consisting  of  two  circular  microstrip  disks  in  a  stacked  configuration 
are  investigated.  Using  a  dyadic  Green’s  formulation,  a  rigorous  analysis  of  the  microstrip 
antenna  is  performed  for  two  stacked  configurations.  Assuming  uniform  current  along  the 
probe,  the  mixed  boundary  value  problem  is  reduced  to  a  set  of  coupled  vector  integral 
equations  using  the  vector  Hankel  transform  and  solved  using  Galerkin’s  method  in  the 
spectral  domain.  Due  to  the  singular  nature  of  the  current  on  the  driven  disk  in  the  vicinity 
of  the  probe,  an  additional  term  is  included  in  the  current  expansion  to  account  for  the 
divergent  nature  of  the  current  near  the  probe  feed  junction  and  insure  continuity  of  the 
current  at  the  junction. 

The  input  imped2mce  eind  radiation  patterns  of  the  stacked  microstrip  antenna  is 
calculated  as  a  function  of  the  layered  substrate,  permittivities  and  thicknesses,  and  the 
ratio  of  the  radii  of  the  two  disks.  Both  dual  frequency  and  wideband  operation  is  discussed. 

Microstrip  discontinuities,  such  as  open  end,  gap  and  step  in  width,  have  been 
widely  studied  by  many  authors.  There  are  different  methods  for  analyzing  microstrip  dis¬ 
continuities,  such  as  quasi-static  approach,  planar  waveguide  model  and  integral  equation 
formulation.  As  the  frequrncy  gets  higher,  the  quasi-static  assumption  is  not  veilid.  In  the 
planar  waveguide  model  analyses,  the  thickness  of  the  substrate  is  assumed  much  smaller 
than  the  wavelength  so  that  a  two-dimensional  model  may  be  applied.  In  this  case,  the 
effect  of  the  radiation  and  the  surface  waves  are  not  considered.  The  integral  equation 
method  has  been  applied  to  study  the  open  end  and  gap  dicontinuities  on  isotropic  sub¬ 
strates.  In  applying  the  integral  equation  method,  various  approximation  were  introduced 
in  the  computation  procedure.  More  recently,  finite  element  expansion  currents  are  used 
to  formulate  a  full-wave  analysis  of  micristrip  discontinuities  on  isotropic  substrates. 


8 


The  open  end,  gap  and  step  in  width  discontinuities  placed  on  anisotropic  substrates 
eire  rigorously  analyzed.  Both  uniaxial  and  tilted  uniaxial  anisotropy  are  considered.  The 
materials  are  assumed  to  be  lossless  and  the  metal  strips  to  be  infinitely  thin.  A  dyadic 
Green’s  function  for  layered  anisotropic  media  is  used  to  formulate  a  set  of  vector  integral 
equations  for  the  current  distribution.  The  fundamental  hybrid  mode  is  assumed  to  be 
propagating  on  the  input  and  output  of  microstrip  lines.  In  solving  the  set  of  vector  integral 
equations,  the  method  of  moment  is  employed.  The  basis  functions  for  the  current  on  the 
metal  strip  consider  the  edge  effect.  Both  logitudinal  ^uld  transverse  currents  are  considered 
in  the  calculation.  The  propagation  constant  for  the  infinitely  long  uniform  microstrip  line 
is  first  calculated.  Then  the  propagation  constant  of  the  fundamental  mode  is  used  to 
formixlate  the  excitation  of  the  discontinuity  problem.  At  the  discontinuity,  loczJ  basis 
functions  are  used  to  simulate  the  local  currents  near  the  discontinuity.  The  scattering 
matrix  can  then  be  obtained,  and  an  equivalent  circuit  model  can  be  proposed.  The  effect 
of  the  anisotropy  is  investigated  and  the  results  are  discussed. 

The  leakage  phenomenon  is  important  in  the  area  of  millimeter-wave  integrated 
circuits  and  integrated  optics.  Theoretical  analyses  and  experiments  have  been  performed 
to  investigate  this  phenomenon.  The  leakage  is  due  to  the  TE-TM  coupling  occurring  at 
the  geometrical  discontinuities,  and  the  leaky  power  in  the  form  of  surface  wave  propagates 
in  the  background  medium. 

There  are  different  methods  to  analyze  the  dielectric  strip  waveguides,  including 
the  approximate  field  matching  method,  effective  dielectric  constant  (EDC)  method,  mode 
matching  method,  etc.  The  first  two  methods  are  approximate,  and  can  not  be  used  to 
predict  the  imaginary  part  of  the  propagation  constant.  In  the  third  one,  ground  planes 
have  to  b**  put  at  some  distance  away  from  the  guiding  structure,  hence  the  effect  of 
radiation  loss  is  neglected. 

An  integral  equation  formulation  using  dyadic  Green’s  function  is  derived  to  solve 
for  the  dispersion  relation  of  single  and  coupled  dielectric  strip  waveguides.  A  method  to 
predict  the  leakage  is  presented,  and  the  leakage  properties  are  investigated. 
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Three  different  dielectric  strip  waveguides  axe  investigated  :  optical  rib  waveguide, 
strip  dielectric  guide,  and  insulated  image  guide.  Both  single  and  coupled  strip  waveguides 
are  studied.  The  cross  section  of  the  dielectric  strips  are  assumed  to  have  rectangular  shape. 
Applying  the  Galerkin’s  method,  the  field  distribution  on  the  cross  section  are  represented 
by  a  set  of  unit  pulse  basis  functions.  Substituting  these  basis  functions  into  the  integrzd 
equations,  and  choosing  the  same  set  of  basis  functions  as  the  testing  functions,  we  cein 
obtain  a  determinant  equation  from  which  the  propagation  constant  can  be  solved. 

For  single  dielectric  strip  waveguide,  it  is  observed  that  the  leakage  occurs  when  the 
effective  refractive  index  is  smaller  than  that  of  a  surface  wave  mode  in  the  background 
medium.  It  is  also  observed  that  if  the  lowest  TE-like  (TM-like)  mode  is  leaky,  the  lowest 
TM-like  (TE-like)  mode  is  non-leaky.  When  the  lowest  order  mode  leaks,  the  surface  wave 
mode  of  opposite  polarization  is  excited.  When  the  higher  order  mode  leaks,  the  surface 
wave  modes  of  both  polarizations  can  be  excited. 

For  two  symmetrical  dielectric  strip  waveguides,  both  the  even  and  odd  modes 
are  investigated.  For  the  leaky  mode,  the  total  leakage  is  due  to  the  leakage  from  each 
individual  strip  waveguide.  At  the  separation  where  the  even  mode  has  a  maximum  leakage, 
it  implies  that  the  surface  wave  modes  excited  by  each  waveguide  add  in  phase.  For  the 
odd  mode  at  about  the  same  separation,  these  coa»al  line  feed,  the  reflection  coefficient 
for  the  TEM  mode  is  obtained  which  allows  one  to  compute  the  input  impedance  at  the 
terminals  of  the  probe.  Numerical  results  for  the  input  impedance  are  presented. 

A  finite  difference  time  domain  technique  for  two  dimensional  time  domain  scat¬ 
tering  of  electromagnetic  waves  is  derived.  The  triangulsir  grids  and  the  control  region 
approximation  are  employed  to  discretize  Maxwell’s  equations.  The  finite  difference  time 
domain  techniques  with  uniform  rectangular  grids  has  been  used  in  the  past.  The  scatter- 
ers  are  modeled  using  staircases  and,  recently,  the  accuracy  of  this  approximation  has  been 
investigated.  Several  types  of  other  grids  have  been  proposed  to  improve  the  staircase  ap¬ 
proximation.  Generalized  nonorthogonal  grid  can  model  scatterer  without  staircasing.  It 
has  been  applied  to  spherical  systems,  yet  they  appear  to  be  cumbersome  for  general  scat- 
terers.  The  “distorted  rectangular  grid”  model  approximates  the  computational  domain 
using  rectangular  grids  and  distorts  the  boundary  grids  to  fit  the  interfaces.  The  triangular 


10 


grid  is  used  in  this  paper,  which  is  very  flexible  in  dealing  with  arbitrary  scatterers  and 
absorbing  boundaries. 

The  control  region  approximation,  which  calls  for  Delaunay  and  Dirichlet  tessella¬ 
tion,  has  been  successfully  applied  to  the  frequency  domain  problems  in  the  past.  Two 
double  integral  terms  are  obtained  by  integrating  the  Helmholtz  equation  about  the  De¬ 
launay  tessellation.  The  lerm  involving  the  Laplace  operator  can  be  converted  to  a  closed 
loop  integral  of  normal  derivatives,  which  can  easily  be  approximated  in  finite  difference 
manner  by  utilizing  the  orthogonal  property  of  Delaunay  and  Dirichlet  tessellation.  The 
remciining  term  can  be  approximated  by  multiplying  the  field  at  the  node  with  the  area.  In 
the  time  domain  problem,  the  same  approximation  is  applied  to  the  wave  equation,  except 
the  term  involving  time  derivatives  is  used  in  time  marching  scheme.  Alternatively,  as  in 
Yee’s  algorithm,  the  first  order  Maxwell’s  equations  are  solved  by  spatially  and  temporally 
separating  the  electric  and  magnetic  fields.  In  the  case  of  electric  polarization,  the  electric 
fields  are  placed  at  the  nodes  and  the  magnetic  fields  are  placed  at  the  center  of  triangular 
edges.  The  curl  H  equation  is  integrated  by  applying  Stoke’s  theorem  and  convert  it  to  a 
closed  loop  integral  of  tangential  magnetic  fields.  This  equation  can  be  used  to  advance 
electric  fields  in  time.  To  update  magnetic  fields,  the  second  curl  equation  is  used.  This 
equation  is  approximated  in  the  finite  difference  manner  by  utilizing  the  orthogonality 
property  of  the  tessellation.  The  equations  for  the  magnetic  polarization  case  can  also  be 
derived  following  the  similar  procedure. 

In  order  to  limit  the  computation  domain,  the  scatterers  are  enclosed  with  artificial 
outer  boundaries.  Continuous  smooth  outer  boundaries,  such  as  circles  and  ellipses,  are 
chosen.  The  second-order  time  domain  absorbing  boundary  conditions  derived  from  the 
pseudo-differential  operator  approach  is  imposed  at  the  outer  boundaries.  These  boundary 
conditions  are  implemented  with  the  control  region  approximation  to  determine  necessary 
field  quantities  at  the  boundary.  The  results  of  the  time  domain  control  region  approach 
are  presented  for  simple  scatterer  geometries,  such  as  conducting  and  coated  cylinders  and 
strips,  by  calculating  both  the  transient  and  time-harmonic  responses. 
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The  Finite-Difference  Time-Domeiin  (FD-TD)  method  was  first  introduced  by  Yee 
who  discretized  Maxwell’s  time  dependent  curl  equations  with  second-order  accurate 
central-difference  approximations  in  both  the  space  and  time  derivatives.  Since  then,  it 
has  been  applied  extensively  to  scattering  and  wave  absorption  problems.  Application  of 
the  FD-TD  method  to  microstrip  problems,  in  which  frequency-domain  approaches  have 
dominated,  has  so  far  attracted  little  attention  until  recently  it  was  used  to  obtsdn  fre¬ 
quency  characteristics  of  microstrip  cavities.  Also,  it  has  been  extended  to  the  analyses 
of  open  microstrip  line  and  microstrip  discontinuity  problems  where  absorbing  boundary 
conditions  are  needed  for  the  simulation  of  the  unbounded  domain.  However,  only  isotropic 
or  simple  anisotropic  media  are  considered  in  the  above  papers. 

A  new  FD-TD  grid  model  is  used  to  solve  microstrip  problems  in  anisotropic  media 
having  tilted  optical  axes  expressed  by  permittivity  or  permeability  tensor  with  off-diagonal 
elements.  This  grid  model  is  indeed  a  superposition  of  two  conventional  grids  with  some 
displacement  which  depends  on  the  optical  axes  of  anisotropy.  Implementations  of  different 
boundary  conditions  are  discussed.  Using  this  model,  the  frequency-dependent  character¬ 
istics  of  microstrip  lines  are  investigated.  The  microstrips  are  assumed  to  be  placed  on 
top  of  anisotropic  substrates  with  tilted  optical  axes.  The  case  with  superstrates  is  also 
investigated. 

In  the  finite  difference  computation,  the  open-end  termination  is  simulated  by  using 
the  open-circuit,  short-circuit  technique.  The  source  plane  is  implemented  by  using  a 
magnetic  wall  with  a  Gaussian  pulse  excited  on  the  surface  under  the  strip.  Because  of  the 
symmetry  of  the  problem,  the  region  under  consideration  can  be  reduced  by  half,  using  a 
magnetic-wall  at  the  center  plane. 

The  fields  at  different  positions  are  first  calulated.  Then  the  Fourier  Transform  is 
taken  to  give  the  field  spectra  from  which  the  voltage  and  current  can  also  be  obtained. 
Using  these  data,  the  effective  permittivity  and  the  characteristic  impedance  can  be  de¬ 
termined.  The  frequency  characteristics  of  microstrip  lines  in  anisotropic  media  obtained 
by  this  method  are  compared  with  the  published  results. 
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Finite  difference  time  domain  (FDTD)  techniques  show  great  promise  in  their  ability 
to  solve  three  dimensional  problems  with  arbitrary  geometry.  Advantages  of  this  method 
include  the  ability  to  model  spatially  or  temporally  varying  media.  These  advantages 
are  due  to  the  complete  discretization  of  both  space  and  time.  Considering  the  volume 
of  information  being  calculated  these  techniques  are  very  efficient  and  are  well  suited  to 
calculation  on  future  parallel  processing  computers.  This  method  was  first  formulated 
by  Yee  in  1966  and  his  basic  algorithm  is  still  in  use.  Recent  work  has  demonstrated  the 
applicability  of  the  FDTD  technique  to  microstrip  problems.  The  centered  finite  difference 
approximations  used  are  second  order  accurate  in  both  space  and  time  yielding  good  results 
for  reasonable  mesh  sizes.  Numerical  techniques  used  to  solve  electromagnetic  problems 
must  limit  the  domain  over  which  the  fields  are  to  be  calculated.  This  mandates  the  use  of 
an  absorbing  boundary  condition  to  simulate  the  outward  propagation  of  waves  incident 
on  the  walls  of  the  mesh.  An  absorbing  boundary  condition  has  been  developed  by  Mur 
based  on  the  work  of  Enquist  and  Majda. 

Our  work  in  this  area  includes  development  of  the  algorithms  mentioned  above  into 
a  general  purpose  computer  code  which  may  be  used  to  solve  for  the  transient  response 
of  electromagnetic  problems  with  an  arbitrary  geometry.  In  addition  to  the  transient 
response,  frequency  domain  parameters  may  be  obtained  by  fourier  transform  of  the  time 
domain  results.  Since  the  fields  are  calculated  throughout  space  and  time  all  other  desired 
parameters  may  be  calc\ilated  from  the  field  quantities.  Specifically,  we  are  analyzing 
rectangular  microstrip  structures  with  as  many  as  two  or  more  ports.  Such  structures  may 
be  used  in  MMIC  filters  or  antennas.  This  problem  is  of  interest  for  several  reasons.  First, 
there  are  existing  frequency  domain  solutions  to  the  resonance  problem  of  a  rectangular 
microstrip  patch,  which  we  may  compare  with  the  FDTD  solution.  Secondly,  the  FDTD 
technique  may  be  used  to  analyze  coupling  of  microstrip  lines  to  the  rectangular  structure. 
This  coupling  may  be  either  a  direct  connection  or  a  gap  coupled  connection.  Advantages 
of  the  FDTD  solution  of  this  problem  are  that  it  is  a  full  wave  solution  which  allows  for 
radiation  or  surface  wave  loss  and  that  no  empirical  values  such  as  ’’effective”  dimensions 
are  needed  for  the  analysis,  also  the  geometry  may  be  altered  easily  to  edlow  for  various 
connections  or  coupling  to  the  patch.  This  is  a  significant  improvement  over  methods 
which  rely  on  a  planar  circuit  approach  in  which  the  substrate  thickness  must  be  smedl 
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compaxed  to  wavelength  and  inherently  three  dimensional  coupling  problems  are  not  easily 
handled.  Comparison  of  our  results  with  various  planar  circuit  approaches  will  be  made. 

A  new  perturbation  series,  coupled  integral  equation  approach  for  calculating  the 
frequency  dependent  circuit  parameters  for  quasi-TEM  transmission  lines  with  lossy  con¬ 
ductors  is  presented.  The  method  considers  the  addition  of  loss  and  dispersion  to  be 
perturbations  on  the  lossless  TEM  case,  and  therefore  the  difference  between  the  propaga¬ 
tion  constant  and  the  wavenumber  in  free  space  is  a  smedl  parameter.  We  obtain  the  lowest 
order  term  of  the  perturbation  series  by  solving  two  quasistatic  problems;  the  electrostatic 
problem  to  get  the  capacitance,  and  the  magnetoquasistatic  problem,  with  the  distribution 
of  current  inside  the  wire  considered,  which  gives  the  frequency-dependent  inductance  and 
resistance.  Both  of  these  problems  are  solved  using  one- dimensional  integral  equations 
for  quantities  on  the  surface  of  the  conductor;  this  represents  a  significant  improvement 
in  efficiency  over  previous  methods.  For  most  cases  of  practical  interest,  the  lowest  order 
term  of  the  series  will  suffice.  If,  however,  the  change  in  the  propagation  constant  from  the 
lossless  case,  due  to  the  altered  inductance  and  the  addition  of  resistance,  is  significant, 
additional  terms  in  the  perturbation  series  can  be  calciolated. 

The  method  is  illustrated  with  the  case  of  one  or  more  wires  embedded  in  a  uniform 
dielectric.  In  the  original  magnetoquasistatic  problem,  the  current  is  entirely  directed  along 
the  cLxis  of  propagation,  and  satisfies  the  frequency- domain  diffusion  equation.  Outside  the 
wire,  the  magnetic  vector  potential  is  in  the  same  direction,  and  obeys  Laplace’s  equation. 
The  boundary  conditions  are  the  continuity  of  tangential  and  normal  magnetic  field  at  the 
interface,  which  can  be  expressed  in  terms  of  the  current  density  and  vector  potential  and 
their  derivatives.  Since  we  can  express  the  ratio  of  the  frequency-dependent  resistance  to 
the  DC  resistance  in  terms  of  the  values  of  the  volume  current  and  its  normal  derivative 
on  the  surface  of  the  wire  only,  we  can  use  a  pair  of  coupled  integral  equations  to  solve  for 
these  quantities  alone,  which  we  can  solve  by  Galerkin’s  method  or  other  finite  element 
methods. 
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Results  obtained  using  this  technique  are  shown  for  some  importeint  cases,  includ¬ 
ing  rectangular  wires,  and  are  compared  with  earlier  methods  and  with  experimental  data. 
Previous  methods  for  calculating  the  resistamce  fall  into  three  categories.  First,  for  cer¬ 
tain  cases,  exact  analytical  results  can  be  obt<uned.  Secondly,  especially  in  the  case  of  a 
rectangular  wire,  the  cross-section  can  be  divided  into  rectangular  segments,  each  much 
smaller  than  a  skin-depth,  across  which  the  current  is  assumed  to  be  constant.  Magne¬ 
toquasistatics  gives  simple  answers  for  the  resistance  and  self-inductance  of  each  element, 
and  the  mutual  inductance  between  elements.  This  leads  to  a  matrix  equation  which  is 
solved  for  the  current  distribution.  The  disadvantages  of  this  technique  are  that  it  requires 
basis  functions  throughout  the  cross-section  of  the  conductor,  which  is  especially  intensive 
as  the  frequency  gets  large.  Also,  closed  form  expressions  for  the  matrix  elements  only 
exist  when  the  elements  are  rectangular  —  other  shapes,  such  as  triangular  patches,  which 
might  be  used  to  fit  a  wire  of  arbitrary  shape,  lead  to  nested  numerical  integrals. 

The  third  method  used  is  a  variational  procedure.  This  is  similar  to  the  method 
presented  here,  except  that  the  current  and  the  magnetic  vector  potential  are  expamded  in 
functions  which  span  the  entire  cross-section.  This  has  two  drawbacks:  first,  it  requires  that 
there  be  a  closed  outer  conductor,  which  is  not  physical  in  many  important  cases.  Second, 
as  in  the  previous  method,  using  elements  which  fill  the  entire  cross-section  increases 
the  computation  time  unnecessarily,  since  only  the  value  of  the  current  and  its  normal 
derivative  at  the  surface  of  the  wire  axe  needed  to  cadculate  the  resistance. 

A  new  method  for  analyzing  frequency-dependent  transmission  line  systems  with 
nonlinear  terminations  is  presented.  The  generalized  scattering  matrix  formulation  is  used 
as  the  foundation  for  the  time  domain  iteration  scheme.  Compared  to  the  admittance 
matrix  approach  proposed  in  a  previous  paper,  it  ha.s  the  advantage  of  shorter  impulse 
response  which  leads  to  sm2Lller  computer  memory  requirement  and  faster  computation 
time.  Examples  of  a  microstrip  line  loaded  with  nonlinear  elements  are  given  to  illustrate 
the  efficiency  of  this  method. 
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As  the  speeds  of  integrated  circuits  increase,  the  effect  of  interconnection  lines  be¬ 
comes  more  and  more  important.  Traditional  lumped  element  circuit  model  must  be 
supplemented  by  the  transmission  line  model  in  order  to  account  for  propagation  de¬ 
lays,  dispersion  and  losses.  This  has  created  needs  for  new  numerical  procedures  that 
can  be  easily  incorporated  into  current  CAD  tools.  To  make  matters  more  complicated, 
the  interconnection  lines  are  terminated  with  not  only  linear  resistors  but  also  nonlinear 
semiconductor  devices,  such  sis  diodes  and  transistors. 

Several  techniques  are  now  commonly  used  to  deal  with  nonlinesu’  circuit  prob¬ 
lems,  for  example,  the  direct  time  domain  approaches,  and  the  semi-frequency-domain  ap¬ 
proaches,  such  as  the  harmonic  balance  technique.  The  semi-frequency-domain  approaches 
are  useful  for  microwave  and  millimeter  wave  integrated  circuits  but  become  impractical 
for  digital  integrated  circuits  because  of  the  latter’s  wide  band  nature.  On  the  other  hand, 
frequency-dependent  psirameters  often  make  it  awkward  to  apply  the  direct  time  domsiin 
approach  to  the  interconnection  line  systems. 

We  propose  a  hybrid  frequency-domain  time-domain  technique  based  on  the  gen¬ 
eralized  scattering  matrix  formulation.  For  an  n  line  system,  we  define  2n  scattering 
parameters  according  to  the  frequency-dependent  characteristic  impedances  of  individual 
lines  {Zoj  =  w)).  The  time-domain  transfer  matrix  (impulse  response)  of  this 

2n-port  system  is  then  obtained  by  the  inverse  Fourier  Transform.  Lastly,  the  nonlin¬ 
ear  equations  associated  with  terminal  characteristics  are  incorporated  and  solved  with 
iteration  procedures  such  as  the  Newton-Ralphson  method. 

The  key  to  efficient  and  stable  solutions  in  this  problem  is  shortening  the  duration 
of  every  transfer  matrix  element.  With  the  generalized  scattering  parameters  approach,  we 
are  able  to  achieve  that  yet  eliminating  the  need  for  artificial  matching  networks  adopted 
by  a  previous  work.  Furthermore,  the  use  of  individual  characteristic  impedances  in  the 
definition  of  scattering  parameters  enables  us  to  gener^llize  this  method  to  coupled  lines 
with  distinct  properties  while  keeping  the  duration  of  transfer  matrix  elements  short.  This 
cannot  be  realized  if  traditional  scattering  parameters  are  used.  We  shall  illustrate  the 
elegance  and  efficiency  of  our  approach  for  a  dispersive  microstrip  line  with  different  non¬ 
linear  loads  and  excited  with  narrow  Gaussian  pulses.  The  elements  of  transfer  matrix 
are  found  out  to  be  either  zero  or  single  retarded  delta-impulse  accompanying  a  small 
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spike  with  very  narrow  spread.  Typical  computation  time  for  a  1000  time-step  iteration 
ranges  from  4  to  27  seconds  on  a  VAXStation  3500.  The  effects  of  line  dispersion  and  load 
nonlinearity  will  be  clearly  delineated  in  the  presentation. 

The  transient  propagation  characteristics  of  VLSI  interconnects  with  discrete  ca¬ 
pacitive  loads  at  various  locations  is  analyzed  based  on  a  hybrid  transmission  lines-lumped 
element  circxiit  model.  Exact  expressions  of  the  Laplace  transform  of  unit  step  responses 
are  first  obtained  through  the  ABCD  matrix  formulation.  We  then  apply  the  equivalent 
dominant  pole  approximation  to  the  transfer  function  with  the  propagation  delays  factored 
out.  The  approximated  transfer  function  can  be  inverted  in  closed  form  and  quickly  eval¬ 
uated.  These  results  provide  efficient  ways  of  finding  approximately  the  effects  on  delays 
and  rise  time  brought  by  VLSI  off-chip  interconnects. 

Because  of  the  dramatic  increase  in  device  densities  on  microelectronic  chips,  the 
propagation  delay  for  off-chip  interconnects  has  become  the  limiting  factor  to  the  speed 
of  VLSI  packages.  Typical  scales  of  these  interconnects  will  be  comparable  or  larger  to 
the  characteristic  wavelength  of  high  frequency  components  of  the  signal.  Therefore,  to 
calculate  the  delays  caused  by  these  interconnects  properly,  a  hybrid  circuit  model  con- 
tsdning  transmission  line  sections  as  well  as  lumped  elements  must  be  used  in  place  of 
the  all-lumped  element  one.  Most  circuit  simulation  packages  are  nevertheless  based  on 
the  latter  and  have  to  resort  to  subsection  approximation  when  deeding  with  transmission 
lines.  This  scheme  will  undoubtedly  lead  to  lengthy  computation,  which  is  not  desirable 
when  a  quick,  heuristic  estimate  of  bounds  are  needed  for  the  initial  phase  of  the  design 
cycle. 

Two  approaches  have  been  developed  for  obtaining  the  approximate  transient  re¬ 
sponse  without  lengthy  simulation.  The  first  kind  of  solution  techniques  emphasize  the 
calculation  of  bounds  to  voltage  responses  from  the  differential  equations  either  by  direct 
integration  or  by  using  the  optimal  control  theory.  On  the  other  hand,  the  second  kind  of 
techniques  analyze  the  properties  of  Laplace  transform  domain  solution.  Thus  far,  their  ap¬ 
plications  are  limited  to  all  lumped-element  and  distributed  RC  networks,  which  can  only 
take  care  of  on-chip  interconnects.  We  shall  take  the  second  approach  by  incorporating 
transmission  line  elements  for  off-chip  delay  estimation. 


17 


Our  configuration  includes  a  series  of  transmission  line  sections  with  arbitrary  dis¬ 
crete  capacitances  and  resistances  loaded  at  junctions.  The  ABCD  matrix  formulation  is 
used  to  obtain  the  Laplace  transform  of  the  unit  step  response.  We  express  the  latter  in 

the  form  of  exp(-5T)/g(«),  where  Q(a)  =  Ao(«)-fi4i(»)exp(*Ti)-f  A3(«)exp(»Tj)-| - with  all 

{Xj(<)}  being  rational  functions  in  a.  The  factor  {exp(-*r)}  is  identified  with  direct  trans¬ 
mission  delay  over  the  total  length  of  the  line.  For  the  rest  part  (1/Q(*)),  we  proceed  to 
apply  the  equivalent  dominant  pole  approximation  technique  [7].  Either  a  single  negative 
real  pole  or  second-order  complex  conjugate  pair  will  be  chosen  for  approximation  depend¬ 
ing  upon  the  property  of  lumped  loads  at  junctions  as  well  as  the  source  impedances.  A 
phase- correction  factor  exp(-«Tm)  is  introduced  to  make  up  for  the  discrepancies  caused 
by  our  low-order  approximation.  The  first-order  and  second  order  approximations  enable 
us  to  obtciin  closed-form  solution  to  the  transient  response.  Comparison  of  the  approxi¬ 
mated  responses  with  those  obtziined  from  brute-force  numericzd  Laplace  inversion  shows 
very  good  match  when  the  propagation  delay  of  an  average  transmission  line  section  is 
less  than  half  the  product  of  junction  load  capacitance  amd  transmission  line  characteristic 
impedance.  Yet  we  only  have  to  spend  a  fraction  of  the  time  for  computations.  The  accu¬ 
racy  of  this  method  will  be  discussed  in  detail  with  some  examples  of  lossless  transmission 
line  networks  in  which  lumped  capacitors  are  loaded  at  regular  intervals. 

The  propagation  properties  of  single  2md  coupled  inhomogeneous  slab  waveguides 
are  analyzed.  An  integral  equation  formulation  using  the  dyadic  Green’s  function  which 
covers  both  the  TE  and  TM  modes  is  proposed.  The  dispersion  relations  are  obtained  by 
applying  the  Galerkin’s  method  to  solve  the  integral  equation.  The  coupling  between  two 
symmetrical  inhomogeneous  slab  waveguides  is  also  investigated.  This  method  is  shown 
to  be  applicable  to  arbitrary  dielectric  constant  profiles. 

The  guidance  and  leakage  properties  of  single  and  coupled  dielectric  strip  waveg¬ 
uides  are  analyzed  using  the  dyadic  Green’s  function  and  integral  equation  formulation. 
Galerkin’s  method  is  used  to  solve  the  integral  equation  for  the  dispersion  relation.  The 
effects  of  the  geometrical  and  the  electrical  parameters  on  the  dispersion  relation  are  inves¬ 
tigated.  A  method  to  predict  the  occurrence  of  leakage  is  proposed.  The  properties  of  the 
even  and  the  odd  leaky  modes  are  also  investigated.  Results  are  compared  with  previous 
analysis  and  shown  to  be  in  good  agreement. 
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A  spectral  domain  dyadic  Green’s  function  for  multilayered  uniajdally  anisotropic 
media  containing  three-dimensional  sources  is  derived.  Tractable  forms  are  shown  to  be 
easily  deduced  from  the  physical  picture  of  the  waves  radiated  by  the  primary  sources 
and  the  multiple  reflections  from  the  stratified  medium.  The  formulation  decomposes 
the  dyadic  Green’s  function  into  TE  and  TM  waves.  The  dyadic  Green’s  function  in  the 
source  region  is  properly  represented  by  extracting  the  delta  function  singularity.  A  simple 
proceedure  to  obtain  the  fields  in  any  arbitrary  layer  is  described.  Recursion  relations  for 
appropriately  defined  reflection  and  transmission  coefficients  are  presented.  Forms  suitable 
for  transmission  line  applications  in  multilayered  media  are  derived. 

Full  mod2J  analysis  is  used  to  study  the  dispersion  characteristics  of  microstrip 
lines  periodically  loaded  with  crossing  strips  in  a  stratified  uniaxially  anisotropic  medium. 
Dyadic  Green’s  functions  in  the  spectral  domain  for  the  multilayered  medium  in  conjunc¬ 
tion  with  the  vector  Fourier  transform  (VFT)  are  used  to  formulate  a  coupled  set  of  vector 
integral  equations  for  the  current  distribution  on  the  signal  line  and  the  crossing  strips. 
Galerkin’s  procedure  is  applied  to  derive  the  eigenvalue  equation  for  the  propagation  con¬ 
stant.  The  effect  of  anisotropy  for  both  open  and  shielded  structures  on  the  stopband 
properties  is  investigated. 

A  direct  three  dimensional  finite  difference  time  domain  (FDTD)  method  is  applied 
to  the  full-wave  analysis  of  various  microstrip  structures.  The  method  is  shown  to  be  ^Ln 
efficient  tool  for  modelling  complicated  microstrip  circuit  components  as  well  as  microstrip 
antennas.  From  the  time  domain  results,  the  input  impedance  of  a  line-fed  rectangular 
patch  antenna  and  the  frequency  dependent  scattering  parameters  of  a  low  pass  filter  and 
a  branch  line  coupler  are  czdculated.  These  circuits  are  fabricated  and  the  measurements 
are  compared  with  the  FDTD  results  and  shown  to  be  in  good  agreement. 
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A  rigorous  dyadic  Green’s  function  formulation  in  the  spectral  domain  is  used  to 
study  the  dispersion  characteristics  of  signal  strip  lines  in  the  presence  of  metallic  crossing 
strips.  A  set  of  coupled  vector  integral  equations  for  the  current  distribution  on  the 
conductors  is  derived.  Galerkin’s  method  is  then  applied  to  derive  the  matrix  eigenvalue 
equation  for  the  propagation  constant.  The  dispersion  properties  of  the  signal  lines  are 
studied  for  both  cases  of  finite  and  infinite  length  crossing  strips.  The  effects  of  the 
structure  dimensions  on  the  passband  and  stopband  characteristics  are  investigated.  For 
crossing  strips  of  finite  length,  the  stopband  is  mainly  affected  by  the  period,  the  crossing 
strip  length,  and  the  separation  between  the  signal  and  the  crossing  strips.  For  crossing 
strips  of  infinite  length  carrying  travelling  waves,  attenuation  cdong  the  signed  line  exists 
over  the  whole  frequency  range  of  operation. 
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Abstract —  In  this  paper,  full  modal  antJysis  is  used  to  study  the  dispersion  character¬ 
istics  of  microstrip  lines  periodically  loaded  with  crossing  strips  in  a  stratified  uniaxialiy 
anisotropic  medium.  Dyadic  Green’s  functions  in  the  spectral  domtdn  for  the  multilayered 
medium  in  conjunction  with  the  vector  Fourier  transform  (VFT)  are  used  to  formulate  a 
coupled  set  of  vector  integral  equations  for  the  current  distribution  on  the  signal  line  and 
the  crossing  strips.  Galerkin’s  procedure  is  applied  to  derive  the  eigenvalue  equation  for 
the  propagation  constant.  The  effect  of  anisotropy  for  both  open  &nd  shielded  structures 
on  the  stopband  properties  is  investigated. 

1.  INTRODUCTION 

In  compact  modules  of  high  performance  computers,  signal  transmission  lines 
between  integrated  circuit  chips  arc  embedded  in  multilayered  dielectric  media. 
These  signal  lines  are  usually  placed  in  different  layers  and  run  perpendicular 
to  each  other.  The  interaction  between  the  orthogonal  crossing  lines  and  the 
signal  bne  affects  its  propagation  characteristics  and  may  cause  considerable  signal 
distortion. 

The  interaction  of  a  pmr  of  crossing  lines  in  an  isotropic  medium  has  been  stud¬ 
ied  using  a  time-domain  approach  [1],  where  coupling  is  described  qualitatively. 
This  method  becomes  computationally  expensive  when  the  number  of  crossing 
bnes  increases.  With  many  identical  crossing  strips  uniformly  distributed  above 
the  signal  line,  the  transmission  properties  arc  characterized  by  stopbands  due  to 
the  periodicity  of  the  structure.  Periodic  structures  have  been  investigated  using 
frequency-domain  methods.  In  [2],  periodically  nonuniform  microstrip  lines  in  an 
enclosure  are  analyzed  on  the  basis  of  a  numerical  field  calculation.  A  technique 
based  on  the  network-analytical  formulism  of  electromagnetic  fields  has  been  used 
to  analyze  striplines  and  finlincs  with  periodic  stubs  [3].  The  propagation  charac¬ 
teristics  of  waves  along  a  periodic  array  of  parallel  signal  lines  in  a  multilayered 
isotropic  structure  in  the  presence  of  a  periodically  perforated  ground  plane  is 
studied  in  [4]  and  that  in  a  mesh-plane  environment  is  studied  in  [5).  More  re¬ 
cently,  the  effect  of  the  geometrical  properties  on  the  propagation  characteristics  of 
strip  lines  with  periodic  crossing  strips  embedded  in  a  shielded  one-layer  isotropic 
medium  have  been  investigated  [6]. 
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In  this  paper,  both  open  and  closed  multilayered  uniaxially  anisotropic  struc¬ 
tures  are  considered.  A  full-wave  aneilysis  is  used  to  study  the  propagation  char¬ 
acteristics  of  a  microstrip  line  in  the  presence  of  crossing  strips.  The  signal  line 
and  the  crossing  strips  are  assumed  to  be  located  in  two  arbitrary  layers  of  a 
stratified  uniaxially  anisotropic  medium.  An  integral  equation  formulation  using 
dyadic  Green’s  functions  in  the  periodically  loaded  structure  is  derived.  Galerkin’s 
method  is  then  used  to  obtain  the  eigenvalue  equation  for  the  propagation  con¬ 
stant.  The  effects  of  anisotropy  on  the  stopband  properties  are  investigated. 
Numerical  results  for  open  and  shielded  three-layer  uniaxially  anisotropic  media 
are  presented. 


II.  FORMULATION  OF  THE  PROBLEM 

In  this  section,  we  present  a  dyadic  Green’s  function  formulation  of  the  problem 
shown  in  Fig.  1(a)  where  the  microstrip  line  and  the  crossing  strips  are  placed 
at  two  different  interfaces  of  a  uniaocially  anisotropic  multilayered  medium.  The 
crossing  strips  are  assumed  to  be  placed  in  a  layer  (t)  and  the  signal  line  to  be 
in  a  layer  (j).  The  crossing  strips  are  considered  to  be  periodic  with  period  p  as 
shown  in  Fig.  1(b).  In  general,  the  permittivity  and  permeability  tensors  of  an 
arbitrary  layer  (1)  are  assumed  to  be  given  by 


where  /  =  0, 1, 2, . . .  ,n, . . . ,  t. 

For  the  stratified  medium,  the  electric  fields  in  layers  (i)  and  {j)  due  to  current 
distributions  Jt(^)  ^.nd  J j(r')  may  be  expressed  as 

£.(f)  =  iw  JJJ^dV  dairy).  Jiif')  +  iu,  jj j^dV' dairy).!  jird)  (3a) 

£j{r)  =  dV  dyry)  •  lii^)  +  dV  f ;,(f,r')  •  7,(r^)  (3a) 

where  r')  is  the  dyadic  Green’s  function  in  layer  (I)  due  to  current  sources 

in  layer  (m). 

For  the  multilayered  structure  shown  in  Fig.  1(a),  the  current  distributions  on 
the  conducting  strips  are  assumed  to  be  surface  currents  lying  on  planes  transverse 
to  z.  Thus,  if  we  are  interested  in  the  transverse  electric  fields  £[*^  and  £^*^  in 
layers  (i)  and  (;),  respectively,  we  can  write 


£<’)(f)  =  ,£;*)(?) + ^!;V) 
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Figure  1.  Geometrical  configuration  of  a  signal  line  in  layer  j  loaded  with 
crossing  strips  in  layer  i. 

where  is  the  transverse  electric  field  in  layer  (1)  due  to  current  sources  in 
layer  (m),  and  is  given  by 

[J^dV’d\’2{fy).Km{r}  (5) 

is  the  surface  current  distribution  in  layer  (m),  and  (2  x  2) 

transverse  part  of  the  dyadic  Green’s  function  G/^(r,r'). 

Since  the  structure  is  assumed  to  be  periodic  in  the  y-direction,  the  electric  field 
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and  the  surface  currents  Km  can  be  expressed  using  the  Floqu-’t  harmonic 
representation  in  the  y-direction.  In  this  case  we  can  get 

_ /  V  coc  rp/2  =(s)  _ 

=  V  (6a) 

J — oo  J  —p/2 


where  I  ~  i,  j,  and 


P  n=-oo-^-“ 

—  2n7r 

r,  =  ix  +  yy,  k,ri  =  +  yPn,  Pn  =  PO  +  - 

P 

Here  p  is  the  period,  Pq  is  the  propagation  constant  of  the  dominant  harmonic 
in  the  Floquet  representation,  and  is  the  spectral  dyadic  Green’s 

function.  ... 

={j)  ={■*)  ={^) 

Using  the  explicit  expressions  for  the  dyadic  Green’s  functions  Gjj  ,  G^j  ,  G^, 

and  Gj/  [Ij,  the  transve  'e  clectiiC  fields  on  the  surface  of  the  conducting  strips 
in  layer  (t)  d":  t-:  ilie  currents  in  layer  (j)  can  be  expressed  in  the  following  form 

(T) 

where  F{ksn,^t)  is  the  kernel  of  the  vector  Fourier  transform  (VFT)  given  by  [8] 

=  A'lc***’'"*  (8) 

and  K j{k,n)  is  the  vector  Fourier  transform  of  the  surface  current  It  is 

given  by 

//I 


The  matrix  ^ij(ksn)  is  given  by 


0  f? 


whose  elements  for  different  i  and  j  are  given  in  Appendix  A. 

In  the  above,  the  transverse  electric  field  expressions  (firn  =  i,j)  satisfy 
the  boundary  conditions  at  the  dielectric  interfaces  of  the  layered  medium.  Apply¬ 
ing  the  final  boundary  condition  that  the  tangential  electric  field  vanishes  on  the 
conducting  strips,  we  can  get  the  following  set  of  dual  vector  integral  equations 
for  the  currents  on  the  metallic  strips 


[  dkx  F{k,„,r,)  ■  ^ii{k,n)  ■  K,{k,n) 

„  =  -oo-'-«’ 

°°  fOO  _  _  _  _  _  _ 

+  Y2  /  dk^  F{k,n,f,)-Uj{k,n)-Kj{k,„)  =  0,  r,  e5.  (11) 

_ J -OO 
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on  the  crossing  strips,  eind 

fOO  _  _  _  _  _  _ 

Y,  /  dk^  F{k,n,r,)-(j,{k,n)-K,(k,r,) 

T.=-oo''-<* 

~/-oo  __ 

+  /  dk^  F{k,n,rs)-ijj{k,n).Kj{k,n)  =  0,  f,eSj  (12) 

n^-ood-°° 

on  the  signal  line. 

The  next  step  is  to  solve  this  coupled  set  of  vector  integral  equations  to  find 
the  dispersion  relation  for  the  signal  line  in  the  presence  of  the  crossing  strips. 

m.  GALERKIN’S  METHOD  AND  THE  EIGENVALUE  EQUATION 

The  formulation  up  to  this  stage  is  exact.  We  now  solve  the  set  of  vector  integral 
equations  (11)  and  (12)  by  using  Galerkin’s  method.  The  unknown  current  distri¬ 
butions  on  the  crossing  strips  Ki{f,)  and  on  the  signal  line  are  expanded 

in  terms  of  the  appropriate  vector  basis  functions  as  follows: 

M  R 

^•mr{x,y)  •  ^mr  (13) 

m=l r=l 

_  K  Q2  _  _ 

Kjiz,y)=Yl  E  (14) 

k=lq=-Qi 

where  Ki{x,y)  and  Kj{z,y)  are  the  surface  currents  on  the  crossing  strips  and 

the  signal  line,  respectively,  $nir  and  are  the  basis  functions,  Amr  and 

arc  the  expansion  coefficients. 

U^g  (13)  and  (14),  the  vector  Fourier  transform  (VFT)  of  the  currents  Ki{r,) 
and  Kj{f,)  are  obtained  as 


M  R 


Fi{k,n)  —  ^  ^ {+)mT,ni^^J0)  '  ■^"»r 

m=l  r=l 

K  Q2  _ 

(15) 

d^j{^sn)  =  ^2  E  ^  (+)kq, ni^x,  P)  •  B  If  g 

k=l  q=-Qt 

(16) 

where 

~  (27r)^  jj  T’’#)  •  ^mr(2:,  y) 

(17) 

and 

-  (2.)2  //  nk.n,Tf.)  ■  'l',(x)e±‘^^v 

(18) 

Substituting  (15)  and  (16)  into  (11)  and  (12),  we  obtain 

M  R  oo  ^  ^  ___  _ 

EE  E  / 

m=:l  r=l  n=— oo 
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;_;1}=-Ql  n=-oo*'“°° 


(19) 


for  rj  €  5{,  and 
M  R  oo 


+ 


M  K  oo  -oo  _  ^  =  _  =  _ 

EE  E  /  ini»'«)  •  ^i^XfP)  ■  •^mr 

m=:l  r=:l  n=— oo 

^  °°  foo  _  _  =_=  _ 

E  E  E  /  F(fc.„,r,).(,.,(A:,„).K(+)  (fcx,/?)-5fc,  =  0 

Jt-l  n=-oo  •'  °° 


(20) 


for  Tj  6  5j. 

=t  .  .  -  _ 

Multiplying  (19)  by  4>„„(s,y)  and  integrating  over  the  support  of  A'j(r,)  for 
u  =  1, 2, . . . ,  M  and  v  —  1, 2, . . . ,  i?,  we  obtain 
M  R  oo 


m  i\.  uu  _  _  _  _ 

m3=l  r=l  n=:— oo  °° 


K  Ql 


'VJ  «J  -oo  ^  =  _  =  _ 

EEE/  » 

t  =  l  g=-Qi  n=-oo  •'  °° 


+ 


1  9=-<?i 


(21) 


Siniilarly,  multiplying  (20)  by  4'j(a:)e~'^‘^  and  integrating  over  the  support  of 
K j{f,)  for  3  =  1,2,. ..  ,K  and  i  =  —Q\,...,0,...,Q2,  we  obtain 

M  R  oo  .oo  ^  ___  _ 

m=lr=ln=— OO  ^ 


i=lo=-D,  n=-oo*'-°° 


+ 


1  <3=-Qi 


(22) 


Equations  (21)  and  (22)  constitute  a  system  of  (5  +  T)  linear  eJgebraic  equations 
with  5  =  MR  and  T  =  K[Qi  +  Q2  +  1)>  ^tid  may  be  written  in  matrix  form  as 


iV-c  =  0 


where 


N  = 


l^nlsxs  [■^i2lsxr 

1^2i1txS  [^22]TxT 


(23a) 

(236) 
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Each  element  ol  the  submatrices  of  N  is  given  by 


°°  t°°  =1  =  _  = 

l^nUmr=  Y.  / 

•/ — oo  '  ’ 

n=-oo  ^ 

(: 

[^12luv,i,=  E 

n——oo  *'  °° 

( 

n=-oo'^~°° 

( 

OO  ^oo 

1^22^5=  E  / 

**—  “OO  *  •* 


For  nontrivial  solution  to  exist,  the  determinant  of  the  coefficient  matrix  of 
(21)  and  (22)  must  be  zero, 


[;v(u.,;3)]  = 


This  is  the  eigenvalue  equation  for  the  propagation  constant  /?  which  describes 
the  dispersion  relation  of  the  loaded  microstrip  line  in  the  multilayered  anisotropic 
medium. 

The  next  step  is  to  choose  appropriate  basis  functions  for  the  surface  currents 
Kj{r,)  and  Ki{r,)  on  the  signal  line  and  the  crossing  strips,  respectively.  The 
expansion  functions  we  use  are 

*  (r  v) _ —  Qr(P)n>i)  0  /2g\ 

<fm,(x,y)-  2^  [  0  Q^(z,Li)J  [  0  ifr(y.u;.)J 


27r  0  Qi,{x,wj) 


where 


^7i(a,7) 

(?n(a,7) 


1  n7r(a  +y/2) 

Rn(0i,y)  =  -  - 

7  7 


(28) 

(29) 

(30) 


w,  and  xvj  are  the  widths  of  the  crossing  strips  and  the  signal  line,  respectively, 
and  L,  is  the  length  of  the  crossing  strips.  When  choosing  the  basis  functions 
for  the  surface  currents,  it  should  be  borne  in  mind  that  the  current  cannot  have 
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a  norm^J  component  to  the  strip  edges.  Futhcrmore,  the  edge  condition  for  the 
parallel  component  should  be  considered.  By  substituting  (26)  and  (27)  into  (17) 
and  (18),  respectively,  we  can  get 

ry,(/3„,u;.)  0  ]  ^  ' 

0  ^(±),(/3n,W,) 


T7  0 


where 


_  4n7r  sin  ^ 

^"("»T')-^(-l)"4-2;2-l-Iy2-  (33) 

Ynia,^)  =  JoK’i  -  l)’^  +  +  ‘^o[(»'  -  l)’^  "  ^l}  (3^) 

=  ^2-2  _-^Xy2  ±  ^7)sm(—  T  — ) 

-  e^(n7r  :^:a7)sin(Y  ±  y)]  (35) 

Equations  (31)  and  (32)  are  then  substituted  into  the  determinantal  equation  (25) 
for  the  calculation  of  the  dispersion  characteristics. 


IV.  NUMERICAL  TREATMENT  AND  RESULTS 

In  this  section,  we  present  numerical  results  for  open  and  closed  three-layer  struc¬ 
tures  with  the  crossing  strips  and  the  signal  line  embedded  in  two  different  layers 
as  shown  in  Figs.  2(a)  and  (b).  In  numerical  cedculation,  the  infinite  series  of 
Floquet  modes  and  the  basis  functions  are  truncated.  The  ranges  of  indices  in 
(24)  are  chosen  as:  —10  <n<9,  A:  =  l,— 1<9<0,  l<m<3,  r  =  l,j  = 
l,-l<t<0,  l<w<3,  and  v  =  1.  It  can  be  seen  that  each  element  in  the  co¬ 
efficient  matrix  can  be  reduced  to  a  sum  of  TE  and  TM  terms,  a  summation  over 
n  Floquet  modes,  and  an  integral  over  kx-  Due  to  the  symmetrical  properties  of 
the  Green’s  function,  the  basis  functions  and  the  test  functions,  all  the  integrands 
are  found  to  be  even  functions  of  kx.  So  the  integration  path  can  be  reduced  to 
an  integral  from  0  to  oo.  In  numeric2J  computation,  the  path  of  integration  in  the 
complex  kx  plane  is  deformed  to  avoid  the  singularities  on  the  real  axis  [9]. 

In  the  following  calculations,  the  parameters  used  for  Fig.  2  are:  di  =  di  =  = 

0.2  mm,  p  =  0.5  mm,  u^]  =  102  =  0.125  mm,  =  1.7  mm,  and  fij  =  =  MO- 

Since  the  crossing  strip  length  is  much  longer  than  the  signal  line  width,  the 
current  near  the  crossing  strip  edges  is  relatively  small  so  that  the  edge  condition 
in  the  basis  functions  can  be  neglected. 
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Figure  2. 


(a) 


(b) 

(a)  Cross  section  of  an  open  structure,  (b)  Cross  section  of  a 
closed  structure. 
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Figure  3  shows  the  dispersion  characteristics  of  a  closed  microstrip  line  loaded 
with  periodic  crossing  strips  (Fig.  2(b)).  The  result  shows  that  the  first  stopband 
appears  due  to  the  coupling  between  the  Floquet  modes  n  =  0  and  n  =  —  1  of 
the  fundamental  mode  of  the  signal  line.  The  upper  and  lower  bounds  of  the 
stopband  is  denoted  by  u>y  and  u>£,  respectively.  At  higher  frequencies,  higher 
order  stopbands  are  encountered  because  of  the  interaction  with  the  higher  order 
modes  of  the  signal  line.  However,  we  concentrate  only  on  the  first  stopband 
which  is  in  the  region  of  practical  interest. 

In  Fig.  4,  the  dispersion  characteristics  of  an  open  and  a  shielded  structure 
are  plotted  in  solid  and  dashed  lines,  respectively.  It  can  be  seen  that  both  the 
stopband  position  and  width  are  close  to  those  of  each  other.  This  is  because  the 
fields  are  mostly  confined  under  the  first  layer  where  the  coupling  between  signal 
line  and  crossing  strips  takes  place.  So  removing  the  top  conducting  plate  does 
not  affect  the  stopband  properties  much  in  this  case.  This  point  is  illustrated  in 
Figs.  5(a)  and  (b)  which  shows  the  effect  of  changing  dj  on  the  stopband  position 
and  width,  respectively.  In  the  following,  we  are  going  to  investigate  the  effect  of 
anisotropy  in  the  second  and  the  third  layer  of  a  closed  structure.  It  is  believed 
that  similar  effects  can  be  observed  in  an  open  structure. 

The  plot  in  Fig.  6  shows  the  effect  of  the  anisotropy  ratio  (AR  =  £2/*2z)  of 
the  second  layer  on  the  stopband  position  and  the  stopband  width.  The  center 
frequency  of  the  stopband  is  not  much  affected  by  the  anisotropy.  However,  the 
stopband  width  is  quite  sensitive  to  it.  The  width  increases  with  l/AR.  For  fixed 
(2,  it  corresponds  to  an  increase  of  «2z,  which  enhances  coupling  between  the 
signal  line  and  the  crossing  strips,  resulting  in  the  rise  of  stopband  width.  For 
1/Ai?  >  1,  it  is  found  in  the  dispersion  diagram  that  a  high  order  stopband  starts 
to  merge  with  the  first  order  stopband,  resulting  in  a  large  stopband  width.  Fig.  6 
is  thus  plotted  up  to  that  value  only. 

In  Fig.  7,  we  investigate  the  effect  of  anisotropy  in  the  third  layer  on  the 
stopband  properties.  As  we  have  expected,  the  stopband  width  is  not  so  sensitive 
to  the  anisotropy  in  the  third  layer  as  it  is  in  the  second  layer  where  coupling 
occurs.  The  change  of  stopband  position  with  the  anisotropy  is  close  to  that  in 
the  second  layer.  Both  are  due  to  the  change  of  the  dispersion  characteristics  of 
the  signal  line  which  results  in  the  lowering  of  the  intersecting  point  of  the  Floquet 
modes  n  =  0  and  n  =  --1.  A  high  order  stopband  is  encountered  for  IjAR  >  1. 

Various  combinations  of  substrate  materials  have  been  used  to  minimize  the 
stopband  width  for  the  closed  structure  (Fig.  2(b)).  The  results  are  summarized 
in  the  following  table: 


Case 

layer  1 

layer  2 

layer  3 

(^0P/’r)c 

A(fcoP/^) 

1 

lOco 

sapph. 

lOeo 

0.3093 

5.93E-3 

2 

lOeo 

10(0 

lOco 

0.3144 

3.67E-3 

3 

■JBHIItM 

O 

o 

0.3061 

3.40E-3 

4 

lOeo 

0.3194 

0.80E-3 

5 

2.3(0 

lOco 

0.3054 

0.53E-3 

6 

«0 

lOeo 

sapph. 

0.3058 

0.87E-3 
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where  (kQpfTr)c  is  the  normalized  center  frequency  of  the  stopband  and  A{kQp/-K) 
is  the  normalized  stopband  width. 

The  two  types  of  anisotropic  substrates  considered  are  Epsilam-10  (e  =  IScq, 
<2  =  10.3eo)  and  sapphire  (e  =  9.4fo,  tz  =  ll.Gco).  Comparison  shows  that  the 
fifth  case  has  the  smallest  stopband  width.  In  fact,  the  stopband  width  is  quite 
sensitive  to  the  separation  of  the  crossing  strips  due  to  the  resonance  effect  [6]. 
Once  the  periodicity  p  is  fixed,  the  stopband  width  can  be  minimized  by  a  proper 
choice  of  substrate  material. 


V.  CONCLUSIONS 

A  dyadic  Green’s  function  formulation  for  the  analysis  of  open  and  closed  mi¬ 
crostrip  lines  in  the  presence  of  periodic  crossing  strips  in  a  stratified  uniaxially 
anisotropic  medium  is  presented.  The  dispersion  characteristics  for  a  three-layer 
structure  is  studied.  Numerical  results  illustrate  the  relationship  between  the 
s*  .ipband  properties  and  the  materied  parameters.  The  effect  of  anisotropy  has 
also  been  investigated.  It  is  found  that  the  crossing  strip  separation  and  the 
anisotropy  in  the  second  layer  are  important  factors  affecting  the  stopband  width. 
To  achieve  small  stopband  width,  careful  choice  of  anisotropy  must  be  made  to 
avoid  the  lowering  of  the  high  order  stopband.  It  should  also  be  noted  that  by  the 
proper  choice  of  substrate  materials,  the  stopband  width  can  be  much  reduced  for 
fixed  crossing  strip  separation. 


APPENDIX  A 

Using  the  dyadic  Green’s  function  formulation  [7],  the  elements  of  ^ij[k,n)  can 
be  obtained  as  follows: 

For  i  =  j,  where  the  source  and  observation  points  are  in  the  same  layer,  we 
have 


,TE  _ 


1 


(i  - 


w 


e  J* 


(Al) 


yTM 

•'JJ 


we. 


‘S:’- 


(A2) 


and  for  j  >  i,  where  the  source  is  in  layer  (j)  and  observation  point  is  in  layer  (i), 
wc  nave 


■TE  _ 


rTE 


(i  - 


(A3) 
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rTM  _  J_ 

»  A:. 


1  - 


vTM  /  \  ,  ,  , 

'i-v  =“!■  "'(1  -  fl?:}')  (i  -  i!3f 


(A4) 

where  and  ®’^®>  respectively,  the  TE  and  TM  upward  transmission 

coefficients  from  layer  (y)  to  layer  (i),  given  by 


vTE  _  vTE  .‘'‘((+1) 


vTM  _  vTM 


I 

1  -  Rif  Mu 


(X5) 


(A6) 


for  I  —  (m  —  2),  (m  —  3), . . . ,  0,  and  for  I  =  (m  —  1),  we  have 


[tTE 

^U(m— l),m 


rTM 


(l  +  rZI) 


(l  +  R^l 

kl 


u(m-l)*  ) 


(l  -  «) 


0- 


u(m-l)  J 


and 


(A7) 


(>18) 


(>19) 


{AW) 
{All) 

Using  the  symmetrical  properties  of  the  dyadic  Green’s  function  in  the  layered 
media,  it  can  be  shown  that  /“^  =  where  a  denotes  TE  or  TM.  In  the  above 
equations,  the  superscripts  {h)  and  (e)  denote  TE  and  TM  fields,  respectively. 

and  R^j  are,  respectively,  the  Fresnel  reflection  coefficients  at  the  lower  and 
upper  boundaries  of  layer  (I)  and  can  be  determined  recursively  by  the  following 
relations 

^{;+i)+:«n(/+i)«  ;  _  ^^^2) 


^  +  -^ru+ij^nu+i) 

where  I  =  0, . . .  ,(n  -  1)  and  R“„  =  R“j. 
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Ki- 


R?n  ,.  +  R°n 


\l-\) 


‘'U(l-l) 


1  +  R?„  ,,R°, 


(>113) 


where  I  =  2, 3, . . . ,  <  and  iijj  =  -Rfo-  are  the  Fresnel  reflection 

coefficients  across  the  interface  between  layers  (/)  and  (/  +  1). 
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Input  Impedance  of  a  Probe-Fed  Stacked 
Circular  Microstrip  Antenna 
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Abstract— The  input  impedance  of  a  microslrip  antenna  consUUnf  of 
two  circular  mkrostrip  disks  in  a  stacked  configuration  driven  b;  a 
coaxial  probe  b  investigated.  A  rigorous  analysis  b  performed  using  a 
dyadic  Green’s  function  formulation  where  the  mix^  boundary  value 
problem  is  reduced  to  a  set  of  coupled  vector  integral  equations  using 
the  vector  Hankel  transform.  Galerkin's  method  b  employed  in  the 
spectral  domain  where  two  sets  of  disk  current  expansions  are  used.  One 
set  b  based  on  the  complete  set  of  ortbogonal  modes  of  tbe  magnetk 
cavity,  and  the  other  employs  Cbebyshev  polynomials  with  the  proper 
edge  condition  for  the  disk  currents.  An  additional  term  is  added  to  the 
dbk  current  expansion  to  properly  model  the  current  in  the  vicinity  of 
the  probe/disk  junction.  Tbe  input  impedance  of  the  stacked  mkrostrip 
antenna  including  the  probe  self-impedance  is  calculated  as  a  function 
of  the  layered  parameters  and  tbe  ratio  of  the  two  disk  radii.  Dbk 
current  distributions  and  radiation  patterns  are  also  presented.  Tbe 
calculated  results  are  compared  with  experimental  data  and  shown  to  be 
in  good  agreement. 


I.  Introduction 

CONVENTIONAL  microstrip  antennas,  consisting  of  a  sin¬ 
gle  conducting  patch  on  a  grounded  dielectric  substrate, 
have  received  much  attention  in  recent  years  [1]  due  to  their 
many  advantages,  including  low  profile,  light  weight,  and  easy 
integration  with  printed  circuits.  However,  due  to  their  resonant 
behavior,  they  radiate  efficiently  only  over  a  narrow  band  of 
frequencies,  with  bandwidths  typically  only  a  few  percent  (!]. 
While  maintaining  the  advantages  of  conventional  single  patch 
microstnp  antennas,  microstrip  antennas  of  stacked  configura¬ 
tions,  consisting  of  one  or  more  conducting  patches  parasitically 
coupled  to  a  driven  patch,  overcome  the  inherent  narrow  band¬ 
width  limitation  by  introducing  additional  resonances  in  the 
frequency  range  of  operation,  achieving  bandwidths  up  to 
10-20%.  In  addition,  stacked  microstrip  configurations  have 
achieved  higher  gains  and  offer  the  possibility  of  dual  frequency 
operation. 

Experimental  work  with  multilayered  microstrip  elements  has 
been  abundant  [2]-I91.  However,  to  date,  theoretical  work  has 
been  relatively  limited,  where  the  study  of  resonant  frequencies, 
modes  and  radiation  patterns  have  been  investigated  (10]-[13]. 
Recently,  the  finite-difference  time-domain  technique  was  ap¬ 
plied  to  stacked  rectangular  microstrip  patch  configurations  [14]. 
There  is  little  or  no  theoretical  analysis  of  the  input  impedance 
of  coaxial  probe-fed  stacked  circular  microstrip  patches.  How¬ 
ever,  the  input  impedance  for  conventional  single-layer  coaxial 
probe-fed  microstrip  antennas  of  circular,  rectangular,  annular 
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Fig.  1.  Stacked  microstrip  antenna  configurations. 


ring,  and  elliptic  geometries  has  been  investigated  by  many 
authors  I15]-I19].  The  impedance  parameters  of  two  planar 
coupled  microstrip  patches  have  also  been  studied  [19],  [20]. 

In  the  calculation  of  the  input  impedance  of  probe-driven 
microstrip  antennas  on  thin  substrates,  the  effect  of  tlie  probe 
results  in  an  additional  inductive  component  to  the  input 
impedance.  This  probe  inductance  has  been  accounted  for  by 
several  authors  through  use  of  a  simple  formula  [19],  [21].  In 
more  rigorous  methods  to  include  the  effects  of  *e  probe,  an 
“attachment  mode”  in  the  disk  current  expansion  is  used  to 
account  for  the  singular  behavior  of  the  disk  current  in  the 
vicinity  of  the  probe,  ensure  continuity  of  the  current  at  the 
probe /disk  junction,  and  speed  up  the  convergence  of  the  solu¬ 
tion.  An  “attachment  mode"  which  represented  the  disk  current 
of  a  lossy  magnetic  cavity  driven  by  a  uniform  cylindrical  probe 
current  was  introduced  in  [16].  More  recently,  similar  and  other 
“attachment  modes,”  with  the  1  /p  dependence  in  the  vicinity  of 
the  probe  and  the  appropriate  boundary  condition  on  normal 
current,  defined  over  the  entire  disk  or  locally  over  a  portion  of 
the  disk,  have  also  been  used  [22] -[25].  In  a  different  approach, 
the  effects  of  the  probe  were  accounted  for  by  expanding  the 
currents  on  the  disk  and  probe  in  terms  of  the  modes  of  a 
cylindrical  magnetic  cavity  satisfying  boundary  conditions  on  the 
eccentrically  located  probe  [26]. 

Considered  here  is  a  microstrip  antenna  consisting  of  two 
circular  microstrip  disks  in  a  stacked  configuration  driven  by  a 
coaxial  probe.  The  two  stacked  configurations  shown  in  Figs. 
1(a)  and  1(b),  denoted  configurations  A  and  B,  respectively,  are 
investigated.  The  disks  are  assumed  to  be  infinitesimally  thin  and 
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perfectly  conducting  and  the  substrates  are  taken  to  be  infinite  in 
extent.  A  rigorous  analysis  of  the  two  stacked  circular  disks  in  a 
layered  medium  is  performed  using  a  dyadic  Green's  fimction 
formulation.  Using  the  vector  Hankel  transform,  the  mixed 
boundary  value  problem  is  reduced  to  a  set  of  coupled  vector 
integral  equations  and  solved  by  employing  Galerkin’s  method  in 
the  spectral  domain.  Two  solutions  using  two  different  basis  sets 
to  expand  the  unknown  disk  currents  are  developed.  The  first  set 
of  basis  functions  used  are  the  complete  set  of  transverse  mag¬ 
netic  (TM)  and  transverse  electric  (TE)  modes  of  a  cylindrical 
cavity  with  magnetic  side  walls.  The  second  set  of  basis  func¬ 
tions  used  employ  Chebyshev  polynomials  and  enforce  the  cur¬ 
rent  edge  condition.  An  additional  term  in  the  current  expansion 
is  taken  to  account  for  the  singular  nature  of  the  current  on  the 
disk  in  the  vicinity  of  the  probe  and  to  ensure  continuity  of 
current  at  the  junction.  This  term,  the  “attachment  mode,”  is 
taken  to  be  the  disk  current  of  a  magnetic  cavity  under  a  uniform 
cylindrical  current  excitation.  It  is  shown  here  explicitly  that 
continuity  of  the  current  at  the  probe /disk  junction  must  be 
enforced  to  rigorously  include  the  probe  self-impedance.  The 
convergence  of  the  results  is  investigated  and  ensured  by  using  a 
proper  number  of  basis  functions.  The  input  impedance  of  the 
stacked  microstrip  antenna  is  calculated  for  different  configura¬ 
tions  of  substrate  parameters  and  disk  radii.  Disk  current  distri¬ 
butions  and  radiation  patterns  are  also  presented.  Finally,  the 
results  are  compared  with  experimental  data  and  shown  to  be  in 
good  agreement.  Throughout  the  analysis,  the  exp(-iwr)  time 
dependence  is  used  and  suppressed. 

n.  Dyadic  Green’s  Function  and  Integral 
Equation  Formulation 

For  a  general  formulation  which  applies  to  both  configurations 
A  and  B  of  Fig.  1,  we  consider  two  coaxial,  circular  perfectly 
conducing  dislu,  of  radii  fl,  and  02,  carrying  current  distribu¬ 
tions  J/f)  =  Kj(^d(z  -  z})  where  j  =  (1, 2)  and  5(-)  is  the 
Dirac  delta  function.  Configuration  A  is  obtained  when  z|  =  </, 
and  Z2  =  d2  and  configuration  B  results  when  zj  =  (Iq  and 
Zi  =  d2. 

Using  the  induced  EMF  method  [27],  a  stationary  formula  for 
the  input  impedance  is  obtained  as 

where  is  the  current  distribution  on  the  probe  and  £  is  the 
total  electric  field  due  to  the  probe  current  and  induced  disk 
currents. 

The  current  on  the  probe,  of  radius  R  and  at  the  position 
Po  ■  (Po-  <^o)’  “  taken  to  be  uniform  and  is  given  by 

d2<z<d2  (2) 

with  local  coordinates  defined  as  Pp  =  p  -  Po"  (Pp>  ^p)- 

Using  a  dyadic  Green’s  function  formulation  in  cylindrical 
coordinates  for  horizontally  stratified  media  [11],  we  obtain 
expressions  for  the  transverse  components  of  the  electric  fields 
due  to  the  disk  and  probe  current  distributions.  Boundary  condi¬ 
tions  require  that  the  transverse  components  of  the  electric  field 
vanish  on  the  perfectly  conducting  disks  and  the  currents  vanish 
off  the  disks,  to  give  the  following  set  of  coupled  integral 


equations  for  the  disk  currents 

[£(p,z  =  z;)]r=  t  e‘"'*rdkpkX{kpP) 

•'0 

=  =^i) 

+  f  rdkpkpj„{k^p) 

m*  — 00  •'0 

•  C.2(*P-  2  =  ^’2)  ■ 

+  f  r dkpkXiKp) 

m--eD  Jo 

=  0,  p<  Oj  (3) 

K'-M-  rdpXikpP) 

Jo 

■^^{kp)  =  0,  p>aj  (4) 

where  j  =  1,2,  and  k^  is  the  transverse  wavenumber  satisfying 
the  dispersion  relation 

kl  +  k}^  =^kj  =  (5) 

in  each  region  /.  In  (3),  z  and  z'  correspond  to  the  longitudinal 
positions  of  observer  and  source,  respectively.  and 

arc  the  vector  Hankel  transforms  of  the  two  disk 
currents  ATilVp)  K^\p),  respectively,  defined  by 

=  r dp pilikpp)  ■  K\i\p)  (6) 

JQ 

where  K^^ip)  is  the  Fourier  coefficient 

Kl^'^{p)  =  ^  rd4>e-‘”%{p)  (7) 

fc  T  J  Q 

and  J„{ki,p)  is  the  kernel  of  the  vector  Hankel  transform 
(VHT)  [28]  given  by 


JUkpP) 

kpP 


-im 

kpP 

J'n,{kpP) 


(8) 


J„{-)  is  the  Bessel  function  of  the  first  kind  of  order  m  and  the 
prime  denotes  differentiation  with  respect  to  the  argument. 

Jm(kpP)  is  the  complex  conjugate  transpose  of  J„(ki,p). 

In  the  last  term  of  (3),  p„(k^)  is  associated  with  the  probe 
current  and  is  given  by 


Pm{kp) 

0 


where 


(9) 


Pn,{kp)-  -:^^J„{k^Po)Jo{kpR)e-‘'”*<^.  (10) 

The  matrix  z)  includes  the  effects  of  the  stratified 

medium  when  relating  the  probe  current  to  the  transverse  elec- 
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trie  fields  and  is  defined  as 


If : 


It  is  clear  that  the  usumed  probe  current  excites  TM  inodes 
only.  The  matrices  z,  zO  with  l,j  =  (1,2)  in  (3)  in¬ 

clude  the  effects  of  the  stratified  medium  when  relating  the  disk 
currents  to  transverse  electric  fields  and  are  of  the  form 


for  m  =  0,±  1,±  2,- •  n=l,2,'",  and  p=I,2,*-- 

correspond  to  the  TM  cavity  modes  and  corre¬ 

spond  to  the  TE  cavity  modes.  The  constants  and  a„p 
correspond  to  the  nth  and  pth  zeros  of  =  0  and 

J^{amp)  =  0,  respectively.  The  VHT  of  these  basis  functions  is 


&kp)  =  0n,Men.n)\ 


0 


{e„Jajf-kl 


32  L.  '^mikpOj) 


The  expressions  for  z),  z,  zO.  and 

ijjikp,  z,  z')  are  given  in  the  Appendix. 

in.  Galerkin’s  Method 

Galerkin’s  method  is  employed  to  solve  the  coupled  vector 
integral  equations  of  (3)  and  (4).  The  currents  on  the  circular 
disks  are  expanded  in  terms  of  a  set  of  basis  functions 

^il>(p)  =  LpKi(p)  +  £6«(p)  (13a) 

n  p 

^2*{p)  =  i:«;(p)  +  i:e*;?i(p)  +  ^j2!«.(p)- 


N  and  P  correspond  to  the  number  basis  functions  V'm>(p)  and 
respectively,  taken  for  the  upper  disk_and  R  and  S 
correspond  to  those  taken  for  the  lower  disk.  /fS!ut(P)  ‘I’c 
“attachment  mode.” 

The  corresponding  VHT  of  the  currents  is  given  by 


-F  (I4a) 

n  p 


A.  TM  and  TE  Modes  of  Cylindrical  Cavities  with 
Magnetic  Side  Walls 

One  set  of  basis  functions  uken  are  those  currents  associated 
with  the  complete  orthogonal  set  of  TM  and  TE  modes  of  a 
cylindrical  cavity  of  radius  aj  (j  =  1,2)  with  magnetic  side 
walls  and  electric  top  and  bottom  walls.  These  current  modes  arc 
given  by 


B.  Chebyshev  Polynomial  Expansion  with  Edge 
Condition 

The  second  set  of  basis  functions  taken  includes  the  edge 
condition  for  the  disk  currents  and  is  taken  to  be  [23] 

P^niploj)  y/l-p^/a^j,  for  p<aj 

0,  for  p  >  fl, 

V  •' 


_  ]iL{p/aj)/ -  P^/°j  >  forp<fl^ 


for  p  >  aj 


for  m  =  0,  ±  1 ,  ±  2,  •  •  •  and  n  =  0, 1 , 2,  •  •  •  T„(x)  is  the 
Chebyshev  polynomial  [29]  and  satisfies  the  recursion  formula 
7'„+i(^)  -  2xr„(x)  +  T„.fx)  =  0  with  Tpix)  =  1  and 

T^(,x)  -  X.  The  term  y/l  -  p^/a^  provides  for  the  proper 
singular  edge  behavior  for  the  azimuthally  directed  current  and 
the  zero  edge  condition  for  the  normally  directed  current.  Since 
the  current  basis  functions  must  have  continuous  current  distri¬ 
butions  on  the  disk,  the  mode  index  m  and  the  Chebyshev 
polynomial  index  n  may  not  be  both  even  or  both  odd  when 
performing  the  current  expansion. 

The  VHT  of  the  above  basis  functions  is  given  by 


ltn\k,)  -  mjt^\k,) 

•fr>Jmn\k^) 


*i/l,{p)  = 


*i^Up)  = 


J'mi^mnPlPj) 


- Ky^mpPlPj) 

npP 

•^mi^mpP  I  Pj) 

0, 


for  p  <  Oj 


for  p  >  Oj 


for  p  <  aj 


for  p  >  fly 


Jt\kp)  -  mjti\k;) 


for  m  &  0.  The  integrals,  with  yj  =  k^aj,  are  defined  by 
(m-n-2\ 

“  8  2  [  \  2  ) 

''^(m  +  n  +  2)/2(yj/2)  •^(m-n-2)/2  (>’>/2) 

I  m  +  n  -  2\ 

~  I  2  j  Am  +  n-2)/2(yj /2) 

■  /[m  -  n  +  2)/2(  .ly  /2) 

+  "'J(m^n)/2 (  yj  „„2  (  >-y  /2)  ( 
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dmn  ~  ^  y  2^ 

(yj/^) 

+  2  d{,m-¥n)P.  (  3'y  /2)  "/(m  -  n)/2  /2) 

~ -A m  +  n  -  2) /2  (  /2  )  ■/(  m  -  n  +  2) /2  {  ^2 /2  )  ] 

(19b) 

■*'"^(Bi+n-2)/2(3y/2)  d{m-n)rii,yjl‘^\ 

(19c) 

dmn  ~  “  7^  [/lm+n)/2(3'2/2)7(„_„)/2(3'2/2)]  • 
yj  ^ 

(I9d) 

In  the  above  expressions,  when  m  and  n  are  not  both  even  or 
both  odd,  the  Bessel  functions  Jm/2(.‘)  are  of  half-integer  order 
130]. 

C.  Attachment  Mode 


D.  Matrix  Equation 

Substituting  the  current  expansion  of  (14)  into  (3),  and  apply¬ 
ing  Parseval’s  theorem,  we  obtain  a  system  of  N  +  P  +  R  +  S 
linear  algebraic  equations  for  each  mode  m  which  may  be 
written  in  matrix  form 

•  Cm  =  d„  (22) 

where 


r 

1  y*\NxN 

r 

1  ^INxP 

1  ^INxR 

[>l]t‘xf 

r  ,*ovn) 

1  -^Ipxn 

r  pOVO) 

1  ■^JPxP 

l-^JPxP 

[>>]$;xf 

r 

1  •^JPxA' 

1  -^JpxP 

[A]ty^ 

[A]tT' 

[A]f:f 

(23) 


The  “attachment  mode”  term  in  the  current  expansion  is 
taken  to  approximate  the  rapidly  varying  currents  in  the  vicinity 
of  the  probe /disk  junction,  ensure  continuity  of  the  current,  and 
speed  up  the  convergence  of  the  solution.  This  term  is  taken  as 
the  disk  current  of  a  magnetic  cavity  of  radius  02  to  a 
uniform  cylindrical  current  source  of  radius  R  positioned  tit  Po 
and  given  by 


T  /n»  -00  J 0  *3^ 


■Jo{l^,R)J„{k,ao) 

ik^I  *  ^im(*-*o) 

4 


im 

T-Mk,p) 

k.p 


Jo{k,R)J„{k,Po)Hi:i^'{k,a^) 

J'mikiaj) 


(20) 

where  //”’(  •)  is  the  Hankel  function  of  the  first  kind  of  order 
m.  The  first  term  in  (20)  is  the  current  induced  on  infinite 
parallel  conducting  planes  by  a  uniform  cylindrical  current.  The 
second  term  is  a  homogeneous  solution_to  the  wave  equation 
added  to  satisfy  the  boundary  condition  H^(p  =  Oj)  =  P*’®" 

viding  for  vanished  normal  current  at  the  edge  of  the  disk. 

The  VHT  of  the  above  attachment  mode  current  distribution 
has  a  closed  form  analytic  expression  given  by 


Uk.p) 

im 


kyP 


Jn,{k^p) 


0  S  p  <  02 


d'mik^aj) 

f  ] 


and 


Each  element  of  the  submatrices  of  A„  is  given  by 

•  lj{k,,z',,z'J)  •  xlilik,) 

Jo 

(26) 


where  yli(^(k^)  and  Xmm(^p)  represent  either  y^l;ii,(k^)  or 
*kmn^k^).  ^ch  element  of  the  excitation  matrix  d„  is  given  by 

dtn  =  -  /  ■  ijy{k^,z',)  ■  Prn{k„) 

Jo 

-  rdk,kyj,]:{k^)  ■  l^(k^,z;,zi)  ■  X«“'(A:J . 
■'0 

(27) 


rv.  Input  Impedance 

Once  the  induced  current  distribution  on  the  microstrip  disks 
due  to  the  coaxial  probe  exciution  is  solved  for,  the  input 
impedance  of  the  stacked  microstrip  antenna  may  be  calculated. 
Applying  (1),  the  input  impedance  for  the  stacked  microstrip 
antenna  is  given  by 

-71/  d<t>pl  dppPp{\E^^{f)]^ 

I  Jdi  •'0  •'0 

+  1^3.2(01*+  (28) 


(21) 
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where  £3/^)  is  die  electric  field  due  to  disk  current  j  and  waveguide.  In  the  small  k^R  limit,  this  term  reduces  to 
is  die  electric  field  due  to  the  probe  current.  After 
integration  over  the  cylindrical  probe  surface  and  some  manipu-  h  J  (k  R)H^'Hk  R) 

lation,  we  arrive  at  *3*^0  4  3  ol  3  0  3 


£  I  dk,k^p*{k^){ri^{k^) 

t  ms  -o»  Jq 

’l3j(^p)[^m  (*p)]p}  +  ^in  whjch  IS  the  formula  used  by  some  authors  as  the  probe  reac- 
2*^  »  <•«  I  tance  [24].  Upon  careful  inspection  of  the  expression  for  the 

=  -  H  /  dk^k^p%{k^)  <  n™(*p)  probe  self-impedance  it  is  noted  that  the  second  term  in 

I  ms  -00  Jq  I  02)  containing  is  zero  when  /?™  is  equal  to  one— the 

N  _  P  case  of  a  probe-fed  parallel  plate  waveguide.  When  is  not 

•  irff'mn[l/'mi(^p)]p  +  ’J3!?(*p)5I*mp  equal  to  one,  this  term  diverges.  This  is  because  the  uniform 

"  current  on  the  probe  leads  to  a  singular  charge  accumulation  at 

P  the  probe  end  giving  rise  to  a  singular  reactance.  Thus,  in  the 

)]  p  case  of  a  microstrip  disk  excited  by  a  probe,  in  order  to  account 

'  properly  for  the  probe-self  impedance,  the  continuity  of  the 

s  _  1  current  at  the  probe-disk  junction  must  be  ensured.  If  we  take 

+>ti^(*p)i:e[«2U*p)]  J  +  the  impedance  due  to  the  probe  and  the  attachment  mode  current 

^  '  ,  ,  together,  we  arrive  at  the  following: 


=  /•— *0^3—  In  {k^R)  =  i60koh^—  In  {k^R) 
2ir  Ho  Mo 


^(selO  ^  2^'“) 


niAk,)  2  k, 

[1  -  !?],“][  1  -  I?7^3][1  -  /?™c'^*3r<‘'o-ri)] 

.  (30) 

tiAKI-  J  [I  -  '  ' 

and  iti  =  y/m/ti-  The  expressions  for  the  generalized  reflec¬ 
tion  coefficients  R^j,  and  /?n/  are  given  in  the  Appendix. 
represents  the  self-impedance  of  the  probe  and  may  be  expressed 

as 

Zllf^  =  ^kM{k,R)H^Ak3R)  - 

4  Jo 


=  ^kM{kyR)Hl'\k,R) 

4 

t  mm-a>  •^m(^2®2) 


where  the  divergent  term  in  the  probe  self-impedance  has 
been  cancelled  by  the  contribution  of  the  attachment  mode  which 
ensures  continuity  of  the  current. 

V.  Radiation  Fields 

The  radiation  field,  or  far  field,  components  in  region  0  may 
be  obtained  from  the  longitudinal  components  with 


fn*  — 


Eon  — 


(32) 


\  *3z  /  For  large  observation  distances,  the  expressions  for  the  field 

components  may  be  evaluated  using  the  saddle  point  method 
Z(2.in)  jj  jjjg  jjjpm  impedance  term  due  to  the  attachment  mode  with  the  saddle  point  being  k^  =  Atq  sin  6  where  6  = 
and  is  given  by  tan"’(p/z).  The  longitudinal  field  components  due  to  disk 


z2"  -  57 

.1  f  , _ /,(*,«)/„(*,>„)  r 

/m--aB  ■'0 


tan"’(p/z).  The  longitudinal  field  components  due  to  disk 
current  A^(p)  arc  given  by 

m«  -oe 


\^3i/  m--» 


« Jeo  tin  *  ^ 


The  first  term  of  the  probe  self-impedance  in  (32)  corresponds 
to  the  input  impedance  of  a  coaxial  probe  driven  parallel-plate 


(36b) 
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where 

/.  ’fi  *1* 

'<■■'(‘'•^'>■777 


[l  - -/?™e'"*i*<M-^i)] 
.gik\^do-z\)g~ikoido 


(37a) 


U  (k  T'\  -  ^ 

#*o  *u 


[l  +  [l  + 

.g-‘kOz<^Ogik\jf.do-t\^ 

V2  * 


(37b) 


(*p.^i)  =  T'? 


2  k2 


[l  -  _  .R™e«*u*i] 

.g-‘k0z‘l0g‘klzktgik2tfi2  (37c) 

*0.2(*P'Zi)  =  "”7^ 

Mo  *2z 
'  [l  - 

.^-M'0j<'0^<*U''1^'*2z*2,  (37d) 

Likewise,  the  radiation  field  component  due  to  the  probe  is  given 
by 


£o''z(p^2)=  L  e'^(-/)’ 


T^^^3(UPp.(*p) 


‘3z 


»i*0' 


kp^kotdit  ^ 


(38) 


where 


and  R  -  0.048<i2-  Given  these  parameters,  the  antenna  is  char¬ 
acterized  by  varying  the  upper  radius  a,  and  the  separation 
between  the  disks  A  2-  Convergent  results  for  the  input  impedance 
and  radiation  fields  using  cavity  mode  basis  fimctions  are  ob¬ 
tained  with  (N  =  4,P  =  3),  {R  =  4,S  =  3),  while  those  using 
the  Chebyshev  polynomial  basis  fimctions  with  the  edge  condi¬ 
tion  use  (N  =  3,P  =  3),  (R  -  3,S  =  3).  For  the  calculation  of 
the  input  impedance  given  by  (29),  very  good  results  are  ob¬ 
tained  with  m  =  ±  1  for  the  terms  associated  with  the  disk 
current  amplitudes,  fiil),,  a^\,  and  and  taking  m  = 
0,±  1,±  2  for  the  probe  self-impedance  and  attachment  mode 
terms  in  (34).  Additional  modes  produce  only  a  1  or  2  Q 
difference  in  the  input  impedance  calculations  for  the  parameters 
considered  here.  It  is  found  numerically  that  the  probe  self-im- 
pedance  and  attachment  mode  impedance  terms  taken  together  in 
(34)  give  rise  to  a  primarily  inductive  reactance  contribution. 
Computation  time  for  the  input  impedance  of  the  stacked  struc¬ 
ture  is  approximately  a  half-hour  of  CPU  time  per  frequency  on 
a  VAXstation  3500. 

Calculated  and  measured  [6]  reflection  coefficients,  F  =  (Z^ 
-Zq)/(Z^  -t-  Zq)  where  Zg  =  SO  0,  are  shown  in  Figs.  2(a) 
and  2(b)  for  the  stacked  configuration  case  A  with  a^  /a2  =  1.01 
and  h2  /a2  equal  to  0.36  and  0.48,  respectively.  The  agreement 
between  the  measured  and  calculated  results  is  very  good.  The 
loop  in  the  impedance  locus  of  Fig.  2(a)  reduces  in  size  in  Fig. 
2(b)  as  />2/a2  increased  from  0.36  to  0.48,  leading  to  the 
wide  bandwidth  behavior  of  this  configuration.  Especially  off 
resonance,  it  is  seen  that  the  attachment  mode  and  probe  self-im¬ 
pedance  terms  are  required  for  accurate  results. 

In  Fig.  3,  return  loss  calculations  using  both  cavity  mode  and 
Chebyshev  polynomial  basis  functions  are  compared  with  mea¬ 
sured  results  [6]  for  the  case  of  Fig.  2(b).  The  agreement 
between  the  calculated  and  measured  results  is  very  good.  A 
frequency  shift  in  the  results  on  the  order  of  1-2%  is  observed 
for  the  two  sets  of  basis  functions.  Due  to  the  presence  of  the 
upper  disk  in  the  stacked  configurations,  two  resonances  associ¬ 
ated  with  the  two  constitutive  resonators  of  the  stacked  structure 
(11]  are  easily  distinguished,  giving  rise  to  a  16  %  —IS  dB 
bandwidth  in  this  case.  One  resonance  is  associated  with  the 
resonator  formed  by  the  lower  disk  and  the  ground  plane  and  the 
second  resonance  is  associated  with  the  resonator  formed  by  the 
two  disks.  Comparing  the  return  loss  of  the  stacked  configura- 


i!  hi  [I  -  *g;i[i  -  «?=)[■  -  «aiii  - 

2  kj  [l  -  _  /?™e'“'**'][l  -  ^ 


(39) 


VI.  Numerical  Results  and  Discussion 

The  integrals  of  the  matrix  elements,  (26)  and  (27),  and  in  the 
impedance  expression  (29),  are  evaluated  numerically  along  an 
integration  path  deformed  below  the  real  axis  to  avoid  the 
singularities  on  the  real  axis  which  correspond  to  the  radiating 
and  guided  modes  of  the  layered  medium.  When  using  cavity 
mode  basis  functions,  the  integrands  vary  asymptotically  as 
1/k^  while  those  using  the  Chebyshev  polynomial  basis  func¬ 
tions  with  the  edge  condition  vary  asymptotically  as  l/k^.  To 
enhance  the  convergence  of  the  integrals  when  using  the  Cheby¬ 
shev  polynomial  basis  functions,  the  asymptotic  values  of  the 
integrands  are  subtracted  out  and  evaluated  analytically. 

For  the  stacked  microstrip  configurations  discussed  here,  the 
following  parameters  are  used:  02  =  1.3233  cm,  h3  =  2  k,  = 
0.115ff2<  <1  *  *3  =  2.45*0,  *2  =  1.22*0  (foam).  Po  ~  0-6ff2. 


tion  to  that  of  the  single  disk,  when  the  upper  disk  and  substrate 
are  removed,  it  is  shown  that  the  input  impedance  of  the  single 
disk  presents  an  impedance  mismatch.  When  the  probe  position 
is  changed  from  0.602  approximately  O.302  to  obtain  a 
match,  a  2.3%  — 10  dB  bandwidth  is  achieved. 

Fig.  4  illustrates  the  effect  of  the  separation  ^2/02  on  the 
input  impedance  for  the  cases  a^|a2  =  105,  where  the  im¬ 
pedance  has  been  calculated  using  the  Chebyshev  polynomial 
basis  functions  including  the  probe  self-impedance  and  attach¬ 
ment  mode  terms.  While  the  position  of  the  lower  resonance 
remains  essentially  the  same,  the  position  of  the  upper  resonance 
is  a  function  of  the  height  decreasing  with  increasing  ^2- 
seen  in  the  figure,  the  excitation  of  the  upper  resonance  in¬ 
creases  with  increasing  h2  (up  to  a  certain  ^2  beyond  which 
there  is  little  coupling  [11]).  This  is  due  to  the  fact  that  the 
coupling  interaction  between  the  two  modes  increases  as  the 


TUUNTSEFF  et.al.:  PROBE-FED  STACKED  CIRCULAR  MICROSTRIP  ANTENNA 


3S7 


0.1  CH2  incrtrn*nt  3.2  -  3.0  CHz 
■  1.01 

O  (3,3)  Ctiebyih«v  (m*t1)  with  ottochmwnt  mod*  (m«0.±1,±2) 

A  (4.3)  Covity  modM  (ma0,±1.±2)  with  attach,  mod*  (ma0,t1,t2) 

•  (3.3)  Chekhov  (matt)  without  ottochmwnt  mod* 

•  Wtoaurwd  [6] 

(a) 


0  t  GHz  incr*m«nt  3.2  -  S.O  GHz 
®i/®j  •  '  01  •  0.46 

o  (3.3)  Ch«byih*v  (matl)  with  ettoehmont  mod*  (ma0.*1.*2) 
o  (4.3)  Cowity  modo*  (ma0,t1,o2)  with  ottoch  mod*  (ma0.t1.*2) 

•  (3,3)  Ch*6y«h*v  (mad)  without  ottochm«rt  mod* 

•  y*o*iwwd  [6] 

(b) 

Fig.  2.  r  of  sucked  configuratioa  A.  a,  /jj  «  l.OI.  (a)  ja-i  «  0,36. 
(b)  *2  /flj  *  0.48. 


upper  resonant  frequency  approaches  the  lower  resonant  fre¬ 
quency,  i.e.,  as  the  upper  resonant  frequency  decreases  with 
increasing  Aj.  Or,  conversely,  the  coupling  interaction  de¬ 
creases  as  the  separation  between  the  disks  approaches  zero. 
Calculated  and  measured  [6]  return  loss  results  are  compared  in 
Fig.  5,  where  a  13%  — 10  dB  bandwidth  is  obtained  in  Fig.  5(a) 
and  where  dual  frequency  operation  is  observed  in  Fig.  5(b). 
Again,  the  agreement  between  the  calculated  and  measured 
results  is  good. 


Frequency  (GHy) 

Fig.  3.  Remm  loss  of  sucked  configuration  A.  a,  ja-^  =■  1.01  and  Aj  jo^ 
B  0.48.  Return  loss  of  single  disk  with  no  upper  substrate  and  Po/^2  ~ 
O.3.O.6. 


Fig.  4.  Input  impedance  of  sucked  cor.fiituralion  A.  0| /Oj  =1.05  and 
*2  /flj  =  0.24  (-),  0.36  (-  -),  0.48  (— ). 


Shown  in  Fig.  6  are  the  input  impedance  results  of  configura¬ 
tions  A  and  B  with  parameters  a,  /Oj  =  1.2  and  Aj -  0.24. 
Generally,  the  two  configurations  have  similar  characteristics. 
For  configuration  B,  the  increased  distance  between  the  two 
disks  and  the  higher  "effective”  dielectric  constant  between  the 
disks  results  in  an  upper  resonance  occurring  at  a  lower  fre¬ 
quency  as  compared  with  configuration  A. 

niustrated  in  Fig.  7  are  the  disk  current  distributions  for 
configuration  A  with  a,  la^  =  1.01  and  Aj/Oj  =  0.48  at  the 
lower  resonance,  that  is  with  AjOj  =  1.655  using  (4,4)  Cheby- 
shev  basis  functions  for  each  disk  and  AjOj  =  1 .68  using  (5,4) 
cavity  mode  basis  functions  for  each  disk.  As  the  number  of 
cavity  mode  basis  functions  is  increased,  the  singular  behavior  at 
the  disk  edge  of  the  i  component  of  the  current  distribution  is 
better  characterized.  The  magnitude  of  the  ^  component  of  the 
current  for  the  upper  disk  is  approximately  uniform  across  the 
disk  where  the  amplitude  slightly  increases  toward  the  edges  due 
to  the  parasitic  effect  of  the  upper  disk  excited  by  the  fringing 
fields. 

In  Fig.  8,  the  radiation  patterns  of  the  stacked  microstrip 
antenna  configuration  of  Fig.  7  are  compared  with  those  of  the 
single  disk  with  no  upper  substrate.  For  the  probe-fed  single 
microstrip  disk,  the  probe  position  is  taken  to  be  Pq/o^  =  0.3, 
while  P0/O2  =  0.6  for  the  stacked  configuration.  TTie  com¬ 
ponent  remains  essentially  the  same  for  both  the  single  disk  and 
stacked  configuration.  The  radiation  pattern  of  the  stacked  con- 
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(a) 


Fig.  7.  Disk  currents  of  stacked  configuration  A.  0,  /Oj  1°  1.01,  /Oj 
O.48.(a)0j|>r^,//l.O»  a^\KjJI\. 


Fig.  5.  Return  loss  of  stacked  configuration  A.  o,  /oj  =  l.OS.  (a) 
=  0.48.  (t)  fij  /02  =  0.36.  (c)  fij  /®j  “  0.24. 


Angle  (degrees) 

Fig.  8.  Radiation  pattern  for  stacked  configuration  A  (— )  with  o,  /Oj  ' 
1.01,  ^2  /®2  “  0.48,  ATjOj  =  1.68,  and  for  single  disk  ( •  •  • )  with 
1.75  and  Po/^2  “  0-3- 


Frequency  (GHz) 

Fig.  6.  Input  impedance  of  stacked  configurations  A  and  B.  0|  /02  »  1.2 
and  *2  /*2  “  0-2^- 

figuration  is  more  directive  than  that  of  the  single  disk,  where 
the  £f  beamwidth  is  decreased  in  the  stacked  case. 

Vn.  Conclusion 

The  input  impedance  of  a  microstrip  antenna  consisting  of  two 
circular  microstrip  disks  in  a  stacked  configuration  driven  by  a 


coaxial  probe  is  investigated.  A  rigorous  analysis  is  performed 
using  a  dyadic  Green’s  function  formulation  where  the  mixed 
boundary  value  problem  is  reduced  to  a  set  of  coupled  vector 
integral  equations  using  the  vector  Hankel  transform.  Galerkin’s 
method  is  employed  in  the  spectral  domain  with  an  additional 
term  used  in  the  current  expansion  to  account  for  the  singular 
nature  of  the  current  in  the  vicinity  of  the  probe,  ensure  continu¬ 
ity  of  the  current,  and  to  speed  up  convergence  of  the  solution. 
Ensuring  continuity  of  the  current  by  means  of  the  attachment 
mode  is  shown  to  be  necessary  for  rigorously  including  the 
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probe  self-inipcdance  and  obtaining  accurate  results  for  the 
stacked  microstrip  configuration.  The  input  impedance  of  the 
stacked  microstrip  antenna  is  calculated  as  a  function  of  the 


lated  results  for  the  stacked  microstrip  configuration  are  shown 
to  compare  well  with  experimental  data. 

Vin.  Appendix 

The  explicit  expressions  for  and 

ijj(.ki,,z,z')  of  (11)  and  (12)  are  given  here.  For  the  stacked 
configurations  of  Fig.  1,  we  have 

[l  -  f?™][l  +  _  e'*3i*3] 

[l  - 

(40) 

[l  - 

’  [l  -  f?^V^*'z'''][l  - 

.g'l‘lz(z\-di)gik2z>‘2  (41) 

where  »)/  =  \/ni/t,  and  zj  =  d2- 
For  source  and  observer  are  in  region  1  (assuming  z  >  z'), 
the  expressions  are 


^  nt  — 


=  0  and  =  - 1  and  f?™  =  1.  The 
coefficients  and  are  defined 
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Abstract-  The  coupled-wave  theory  is  generalized  to  analyze  the  diffraction  of  waves 
by  chiral  gratings  for  arbitrary  angles  of  incidence  and  polarizations.  NumerierJ  results 
are  illustrated  for  the  Stokes  parameters  of  diffracted  Floquet  modes  versus  the  thick¬ 
ness  of  chiral  gratings  with  various  chiralities.  Both  horizontrd  and  vertical  incidences 
are  considered  for  illustration.  The  diffracted  waves  from  chiral  gratings  are  in  general 
elliptically  polarized;  and  in  some  particular  instances,  it  is  possible  for  chiral  gratings 
to  convert  a  linearly  polarized  incident  field  into  two  nearly  circularly  polarized  Floquet 
modes  propagating  in  different  directions. 

1.  INTRODUCTION 

There  has  been  considerable  interest  in  the  theoretical  study  of  scattering  from 
chiral  media.  A  detailed  account  regarding  the  chiral  medium  (or  optical  activity) 
and  some  suggested  applications  can  be  found  in  [1-3].  Chiral  medium  charac¬ 
terized  by  a  biisotropic  constitutive  relation  is  a  special  case  of  the  bianisotropic 
medium  whose  electromagnetic  properties  have  been  extensively  studied  by  Kong 
[4,5].  Jaggard  et  al.  [6]  studied  the  propagation  of  electromagnetic  waves  through 
a  random  collection  of  short  helices  and  demonstrated  the  physical  basis  of  chiral¬ 
ity.  Recently,  Bassiri  et  al.  [1]  analyzed  the  reflection  and  transmission  of  waves 
from  a  half-space  chiral  medium.  Viitanen  et  al.  [7]  solved  the  eigensolutions  for 
the  reflection  of  waves  by  the  interface  of  two  chiral  half-spaces.  Jaggard  et  al.  |2] 
analyzed  the  propagation  of  electromagnetic  waves  in  periodic  chir0il  structures 
under  the  weak  coupling  approximation  for  the  case  that  the  wave  is  normally 
incident  on  a  chiral  slab  with  the  fringes  parallel  to  the  interface. 

Periodic  gratings  have  also  been  the  object  of  extensive  research  through  the 
years  because  of  its  many  applications  in  distributed  feedback  laser,  integrated 
optics,  acousto-optics,  quantum  electronics,  and  holography.  For  the  analysis  of 
wave  diffraction  by  periodic  surface  gratings,  methods  including  the  method  of 
moments  [8,9]  and  extended  boundary  condition  method  [10]  are  rigorous  and  in 
general  computationally  efficient.  For  the  analysis  of  periodic  slanted  dielectric 
gratings,  a  rigorous  coupled-wave  method  has  been  developed  [11]  and  generalized 
to  anisotropic  gratings  by  [12,13]. 

In  this  paper,  we  present  a  generalization  of  the  coupled-wave  method  to  inves¬ 
tigate  the  diffraction  of  electromagnetic  waves  by  a  periodic  chiral  grating  with 
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the  fringes  allowed  to  be  slanted.  The  coupled-wave  approach  writes  the  fields 
inside  the  grating  region  by  Fourier  scries  expansions  and  converts  the  Maxwell’s 
equations  into  a  set  of  coupled  first-order  differential  equations  for  the  Fourier  co¬ 
efficients.  This  system  of  equations  is  solved  by  the  state  variable  method  which 
gives  the  eigenmodes  of  fields  inside  the  grating.  Subsequently,  the  fields  inside 
the  grating  are  expanded  by  these  eigenmodes  and  the  unknown  coefficients  of  the 
diffracted  Floquet  modes  outside  the  grating  are  solved  by  matching  the  boundary 
conditions  on  the  interfaces. 


11.  COUPLED  WAVE  THEORY 


Consider  the  configuration  (Fig.  1)  that  a  plane  wave  with  the  time  harmonic 
dependence  [exp(— impinges  on  a  periodic  chirad  grating.  The  electric  field 
of  the  incident  wave  is  given  by 


Ei  =  iicxp{ki  ■  r) 

having  the  incident  wave  vector  Jbj  defined  as 

fcj  ^  “h  y^yi 
k^i  =5  k  sin  6  cos 
kyi  =:  k  sin  6  sin 
kzi  =t  k  cos  0 

where  k  is  the  wavenumber  in  free  space. 

Z 


(1) 


(2) 

(за) 

(зб) 
{3c) 


i 


f 


Figure  1.  Configuration  of  the  diffraction  from  a  slanted  modulated  chiral 
grating. 


and  expanding  the  diffracted  fields  outside  of  the  grating  region  by  the  Floquet 
modes  yield  the  following  expressions  for  the  total  fields  above  and  below  the 
grating  region,  E  and  E2,  respectively. 


E  =  Ei  +  Y^RnCXp^ikn-r),  2>0 
n 

(5a) 

E2  =  ^  Tn  exp  {-inGzd)  exp  |tfc2n  •  +  dz)]  , 

z  <  —d 

(56) 

where 

kn  =  xkxn  +  yk„i  +  zkxn 

(6) 

k^n  “  kkxn  “b  ifkyi  zfcjxn 

(7) 

and  kxn,  ktn  ,  and  k2zn  are  given  by 

kxn  —  kxi  "b  riGi 
kzn  =  yjk^  —  k\,^  -  k^ 

k2zn= 

(8a) 

(86) 

(8c) 

where  the  branch  cut  of  the  square  root  operation  is  chosen  along  the  negative 
real  axis  of  its  argument. 

Assuming  that  the  medium  in  the  grating  region  is  characterized  by  the  con¬ 
stitutive  relations  for  chiral  medium  [14] 

D  =  fQtE (9u) 
B  —  H  -  ixcy/m^o  E  (9^) 

where  «  and  n  are  the  relative  permittivity  and  permeability,  and  Xc  charac¬ 
terizes  the  chirality  or  handedness  of  the  chiral  medium.  Xc  is  a  dimensionless 
quantity  and  Xc  1  denotes  smEill  chirality.  The  constitutive  relations  given  by 
(9)  can  be  related  to  the  constitutive  relations  [2]  by  the  following  transformations 

Xc  ^cM\/mo/<0 

where  fc  is  the  chirality  admittance. 

In  the  constitutive  relations  given  by  (9),  t,  Mi  ®-nd  are  all  assumed  to  be 
periodic  in  the  direction  of  grating  vector  G  so  that  the  following  Fourier  series 
expansion  of  the  grating  parameters  can  be  defined  as 


=  y~]<n  exp  (mg -f) 

(10a) 

=  ^  Mn  exp  {inG  ■  r) 

(106) 

=  ^  Xn  exp  {inG  •  f ) 

(10c) 
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■  '■) 

Xc  -  ^ 

-Y^ -  =  13 

Xc  -  M«  ^ 


X?  - 


^Cn  exp(tnGT) 


Substituting  the  constitutive  relations,  (9a)  and  {9b),  into  Maxwell’s  equations 
yields 

V  X  £  =  iknTjU  +  kxcE  (11) 

V  X  rjH  =  —iktE  +  kxc'nH  (12) 

where  q  =  represents  the  free  space  impedance. 

For  the  coupled  wave  theory,  electromagnetic  fields  inside  the  grating  region 
(0  >  2  >  —  d)  are  assumed  to  have  the  following  series  expansions 

£  =  13  [•^in(2)  X  +  Eyn{x)  y  +  Ezn{x)  i]  exp  (t/?„  •  r)  (13a) 

£  =  13  [Hzn{z)  X  +  Hyn{z)  y  +  H,n{x)  x]  exp  (t/?„  •  r)  (136) 

where 

/?„  =  ikin  +  ykyi  +  znGz  (14) 

Substituting  (10)  and  (13)  into  (11)  and  (12)  yields  the  following  coupled  wave 
equations 

^ 

~  zEzn  +  ik  ^  ^  y-n—mTlHym  +  k  ^  ^  Xn—mEym  “i"  ^k^nExn  (15a) 

m  m 

I ^  =  ~  inGxEyn  “  ik  “  k  Xn~mExm  +  iky^Ezn  (15^>) 


—  —  iTiGzV^m  “  ^k  ^  ^  ^n^mEym  "I"  k  ^  Xn^-^mV^ym  “i"  ^kxn^E zn  (15c) 


drjHyn 


_  —  inGz'n^yn  +  ^k  ^n—mEzm  ~~  k  ^  ^  Xn—m^^xm  4"  ^ky^TjHzn  (15rf) 


Ezn  —  ^  ^  *^n— 


Eym  ~  "7 —  Exm  )  Tn— 


iHzn  =  13  -<n-m  ~  ^  ^  q£*m) 


In  deriving  (16),  the  z  components  of  E  and  H  are  first  solved  in  terms  of  the 
z  components  of  V  x  £  and  V  x  H  using  (11)  and  (12).  Subsequently,  (16)  is 
obtained  by  substituting  (13)  into  the  resulting  equations.  Note  that  the  chirality 
of  the  gratings  can  be  considered  small  in  (15a) -(15d)  if 

«.  M  ~  0(1)  »  Xc 
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On  the  other  hand,  if  Xc  is  of  the  order  of  one,  then  it  would  have  significant 
amount  of  contribution. 

By  substituting  (16)  into  (15),  (15a) -(ISd)  can  be  recast  into  a  matrix  form 


as 


where 


—  V  =  AV 
dz 


r 


V  = 


■>x/n 


f)Hxn 

L’l-^yn 


(17) 


(18) 


and  matrix  A  is  the  coupling  matrix.  Note  that  in  the  above  equation  vector  V 
shovild  be  realized  as  containing  the  modes  with  all  indices  n’s  and  in  fact  is  an 
infinite  dimensional  vector. 

The  solutions  to  the  coupled-wave  equations  can  be  expressed  as 

["■f’l.nml 


V  —  ^  ^  Cmf^m  exp  (.^mx)  —  ^  ^  C-m 


P'1,  nm 

Pi,  nm 


exp(Am2) 


(19) 


P4,nm 


where  Pm  is  the  eigenvector  corresponding  to  the  eigenvalue  Am  for  the  matrix 
Pl,nmi  P2,nmt  P3,nm  i  and  P^^nm  correspond  to  £xn,  Pyny  V^znt  and  T)Hyn, 
respectively.  In  practice,  the  series  expansions  will  be  truncated  so  that  matrix 
A  becomes  a  finite  matrix.  Eigenvalues  and  eigenvectors  of  the  truncated  ma¬ 
trix  are  solved  by  using  a  computer  library  program  EIGCC  in  the  International 
Mathematics  and  Statistics  Library  (IMSL).  Cm’s  are  the  unknown  coefficients  to 
be  determined  from  boundary  conditions. 

Note  that  Gauss’  law  for  the  source-free  region  gives  two  additional  equations 

*  Pn  —  0  (20a) 

kin  ■Tn  =  0  (206) 


Furthermore,  the  continuity  of  ihe  tangential  components  of  electric  and  magnetic 
fields  on  the  upper  and  lower  interfaces  of  gratings  yields  the  following  equations: 


15*n(0)  = 

+  Pzn 

(21) 

o 

II 

+  Pjfn 

(22) 

krjHxniO)  = 

Snoif^yi^z  +  ^liCy)  + 

{^■yiRzn  ~  ^znRyn) 

(23) 

kvHyniO)  = 

6no(~^x«®i  —  ^zi^z) 

4*  (^zn^zn  “  ^zn^zn) 

(24) 

•£'in(“d)  = 

T^n 

(25) 

1 

II 

Tyn 

(26) 

kr^H^ni-d)  = 

kyiTxn  +  kiznPyn 

(27) 
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—  ~^2zn^xn  ~  ^xn^zn 


(28) 


where  6„m  is  the  Kronecker  delta  and  e*,  Cy ,  and  e^  are  the  x,  y,  and  z  com¬ 
ponents  of  Cj  which  defines  the  polarization  of  the  incident  electric  field  given  by 
(!)• 

The  above  equations  (20)-(28)  can  be  used  to  solve  all  the  unknowns.  For 
example,  if  N  modes  are  kept,  there  are  47V  eigenvectors  inside  the  grating 
region  and  37V  components  for  each  of  Rn  and  Tn-  Note  that  the  total  number  of 
unknowns  is  ION  which  equals  the  total  number  of  equations.  For  the  efficiency 
of  numerical  computation,  the  components  of  Rn  and  Tn  arc  first  expressed 
analytically  in  terms  of  Cm’s.  Thus  (20)-(28)  reduce  to  47V  equations  for  c„i’s. 

After  Tn  is  solved,  the  horizontal  and  vertical  components  of  the  nth  trans¬ 
mitted  Floquet  mode  are  defined  as 


Eh  =  h-Tn  (29a) 

E^  =  v-Tn  {29b) 


where 


-  -kyiX  +  kxnp 
- - 

•  ^2in  kxn^  +  ky^y  kp  . 
k2  kp  ^  k2 


(30a) 

(30b) 

(31) 


and  1:2  ^he  wavenumber  of  lower  half-space. 

Similarly,  the  horizontal  and  vertical  components  of  the  nth  reflected  Floquet 
mode  are  defined  as 


Eff  —  h  •  Rn 
Ev  —  V  •  Rn 

where  k  is  still  defined  by  (30a)  and  v  is  given  as 

-  kzn  kxnk  -f-  kyiy  kp 

V  = - - —  -) - -1 

I  T  '  r  ' 


(32a) 

(326) 

(33) 


The  diffraction  efficiency  is  defined  as  the  ratio  of  the  z  component  of  the 
power  carried  by  Floquet  modes  to  the  z  component  of  the  power  of  the  incident 
wave. 

i2 


D Ern  —  Rc  (^rn/l:zt)|7Zn| 


(34a) 


DEin  =  R<^{k2xn/kzi)\Tn\^  (346) 

For  lossless  gratings,  the  sum  of  the  efficiencies  for  aU  the  propagating  modes 
should  be  unity.  For  the  numerical  examples  illustrated  in  this  paper,  the  coupled 
wave  method  gives  highly  convergent  results  where  the  computed  total  diffrac¬ 
tion  efficiency  deviates  from  1  by  0.02%  for  the  worst  case.  Note  that  a  VAX 
workstation  is  used  for  our  numerical  computation. 
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III.  NUMERICAL  RESULTS  AND  DISCUSSIONS 

In  this  section,  numerical  examples  for  the  effects  of  chirality  are  illustrated  for  the 
polarization  slates  of  the  waves  diffracted  by  gratings  with  rectangular  grooves 
(Fig.  2),  even  though  the  theory  developed  in  Section  II  can  be  applied  for  much 
more  general  cases.  However,  gratings  with  sinusoidal  modulation  have  also  been 
anedyzed  but  the  features  observed  are  essentially  the  same  as  those  of  rectangular 
modulations  presented  in  this  paper. 


Z 


Figure  Z.  Wave  diffraction  from  a  rectangular  chiral  grating. 

Figure  2  shows  the  geometry  of  the  rectangular  grating  for  the  cases  presented. 
For  all  the  cases  considered,  P  equzJs  A  which  is  the  electromagnetic  wavelength 
in  free  space.  The  relative  permittivity  of  the  substrate  below  the  corrugated 
region  is  «2  =  2.5,  and  the  relative  permeability  is  one  for  all  regions.  The 
chirality  parameter  (xc)  is  chosen  as  0,  0.188,  and  0.377  to  illustrate  the  effects 
of  chirality. 

Two  sets  of  incident  angles  (0,  <f>)  are  selected:  (30  deg,  0  deg)  and  (35.26  deg, 
30  deg).  The  former  one  corresponds  to  the  situation  that  the  incident  direction 
is  perpendicular  to  the  ruling  direction;  wherezis  the  latter  one  is  not.  In  either 
case  the  incident  angle  is  chosen  such  that 

kxn\  —  ~kxi  (35) 

ln=— 1 

In  other  words,  the  wave  is  incident  at  the  first  Bragg  angle.  For  the  cases  studied 
in  this  paper,  most  of  the  power  are  carried  by  the  —1  and  0  transmitted  Floquet 
modes. 
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Hence  only  the  diffraction  efficiency  and  polarization  state  of  these  two  domi¬ 
nant  modes  are  illustrated.  In  order  to  present  the  polarization  state  of  diffracted 
Floquet  modes  in  a  compact  format,  orientation  (o)  and  ellipticity  (/3)  angles 
(Fig.  3)  arc  utilized  to  express  the  transmitting  and  receiving  polarization  states. 
Here,  the  definitions  from  (l)  are  adopted.  Therefore,  horizontal  (H)  and  vertical 
(V)  polarization  states  will  have  zero  degree  ellipticity  angles,  with  orientation 
angles  of  Odeg  and  90 deg.  Right  and  left  polarization  states  are  orientation 
independent  with  ellipticity  angles  of  45  deg  and  —45  deg,  respectively.  In  addi¬ 
tion,  Odeg  <  Q  <  180 deg  and  —45 deg  <0  <  45 deg. 

POLARIZATION  ELLIPSE 


V 


Figure  3.  Polarization  ellipse,  a  and  0  represent  orientation  and  ellip¬ 
ticity  angles,  respectively. 


A  general  polarization  vector, 
P  = 


\Eh\txp{i4>hy 

_\Ev\^x'p{i<Pv) . 

(36) 


written  in  terms  of  horizontal  (h)  and  vertical  (r)  vector  components,  can  be 
transformed  into  a  normalized  Stokes  vector  as  follows 


(37) 


where  Using  this  equation,  the  angles  a  and  0  are  obtained. 


•50  ■ 

■  ' 

1 

51 

52 

= 

-  \Ev\^ 
2|£/,||£B|cosd' 

=  (|£^P  +  |£.|2) 

cos  20  cos  2a 
cos  20  sin  2a 

.53. 

.  sin  (i>  _ 

sin  20 
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Figure  4.  Diffraction  efficiencies  and  Stokes  parameters  of  transmitted 
Floquet  modes  0  and  —1  for  horizonteil  polarization  incidence. 

0  =  30  deg  and  4>  =  Odeg.  (a)-(c):  xc  =  0;  (d)-(f):  Xc  = 
0.188;  (g)-(i)'.  Xc  =  0,377.  Solid  and  dashed  lines  correspond 
to  Floquet  mode  0  and  —1,  respectively. 

Figures  4  and  5  illustrate  the  Stokes  parameters  of  the  transmitted  0  and  —  1 
Floquet  modes  versus  the  thickness  of  rectangular  gratings  (d),  respectively,  for 
horizontally  and  vertically  polarized  incidences.  The  incident  angle  is  0  =  30  deg 
and  4>  =  Odeg.  For  the  case  of  zero  chirzility  (achiral  grating)  and  horizontal 
incidence  [Fig.  4(a)-(c)l  which  has  been  studied  in  [15],  the  diffraction  efficiencies 
oscillate  when  the  thickness  of  the  grating  increases.  Note  that  for  Xc  =  0  the 
diffracted  polarizations  remain  horizontally  polarized  (q  =  Odeg  and  P  =  Odeg) 
in  Figs.  4(a)-(c)  and  vertically  polarized  (a  =  Odeg  and  P  =  Odeg)  in  Figs. 
5(a)-(c).  When  the  grating  material  becomes  chiral  [xc  0)i  ^he  diffracted 
Floquet  modes  are  in  general  elliptically  polarized. 
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Figure  5.  Diffraction  efficiencies  and  Stokes  parameters  of  transmitted 
Floquet  modes  0  and  —1  for  vertical  polarization  incidence. 

0  =  30  deg  and  4>  =  0  deg .  (a)-(c):  Xc  =  0;  (d)-(f):  Xc  = 
0.188;  (g)-(i):  Xc  =  0.377.  Solid  and  daished  lines  correspond  to 
Floquet  mode  0  and  —1,  respectively. 

Figures  6  and  7  represent  the  cases  that  the  incident  direction  is  not  perpen¬ 
dicular  to  the  ruling  direction  {4>  =  30  deg  and  0  =  35.26  deg),  respectively,  for 
horizontal  and  vertical  incidence.  All  the  diffracted  modes  arc  generally  ellipti- 
cally  polarized.  In  contrast  to  the  cases  of  ^  =  Odeg  [Figs.  4  and  5  (a)-(c)].  Figs. 
6(c)  and  7(c)  show  that  the  Floquet  modes  are  ellipticaUy  polarized  for  achiral 
gratings.  This  is  due  to  the  fact  that  horizontal  and  vertical  polarizations  are 
coupled  for  an  oblique  incidence  (i^  ^  0). 

Several  features  regarding  the  effects  of  chirality  can  be  observed  in  all  cases. 
The  chirality  introduces  more  coupling  between  0  and  —1  modes,  thus  resulting  in 


Mitd- 
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more  even  distribution  of  the  power  between  these  two  modes,  therefore,  smaller 
fluctuation  for  the  diffraction  efficiencies  as  a  function  of  thickness  d.  In  addition, 
like  the  diffraction  efficiencies,  the  ellipticity  angles  also  show  oscillations  [(f)  and 
(i)  of  Figs.  4-7]  as  the  chirality  is  introduced.  Moreover,  at  some  particular 
instances,  it  is  possible  for  chiral  gratings  to  convert  a  linearly  polarized  incident 
field  into  two  nearly  circularly  polarized  Floquet  modes  propagating  in  different 
directions.  For  example,  for  the  case  of  Xc  =  0.188  ((c  =  0.0005  mho),  when  the 
depth  of  the  grating  is  about  2.4  A,  both  modes  are  nearly  circularly  polarized 
(0  =  ±45  deg). 


(b)  a:  H  ioc.,  Xt  *0  377 


(0  fi:  H  iDC..  Xo  -  0.377 


Figure  6.  Diffraction  efficiencies  and  Stokes  parameters  of  transmitted 
Floquet  modes  0  and  —1  for  horizontal  polarization  incidence. 
B  =  35.26 deg  and  4>  =  30 deg.  (a)-(c);  Xc  =  0;  (d)-(f): 
Xc  =  0.188;  (g)-(i):  Xc  =  0.377.  Solid  and  dashed  lines  cor¬ 
respond  to  Floquet  mode  0  and  —1,  respectively. 
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Figure  7.  Diffraction  efficiencies  and  Stokes  parameters  of  transmitted 
Floquet  modes  0  and  —1  for  vertical  polarization  incidence. 

6  —  35.26 deg  and  4>  =  30 deg.  (a)-(c):  Xc  =  0;  (d)-(f): 
Xc  =  0.188;  (g)-(i):  Xc  =  0.377.  Solid  and  dashed  lines  cor¬ 
respond  to  Floquet  mode  0  and  —1,  respectively. 

Note  that  the  orientation  angle  of  the  diffracted  waves  in  general  has  a  nef 
ative  slope  versus  thickness  d  of  the  grating  in  the  cases  studied  {xc  >  0  c 
right-hand  chiral).  In  fact  when  chirality  Xc  becomes  negative  (left-hand  chiral 
the  curves  of  orientation  angles  become  the  mirror  images  of  the  correspondii 
right-hand  chiral  cases  with  respective  to  o  =  90  deg  .  The  phenomenon  observe 
can  be  explained  by  studying  how  the  orientation  of  waves  rotates  as  it  propagat 
through  a  chiral  slab  (3).  Neglecting  the  reflection  from  the  boundaries,  a  lineai 
polarized  wave  impinging  on  a  chiral  slab  can  be  divided  into  two  circularly  pol 
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ized  waves  of  opposite  handedness.  For  the  case  of  right-hand  chiral  medium,  the 
right-hand-circularly  polarized  wave  travels  slower  than  the  left-hand-circularly 
polarized  wave.  If  we  trace  the  polarization  of  the  total  field  along  the  propa¬ 
gation  path  at  any  given  time,  then  the  resulting  polarization  trajectory  would 
circulate  in  the  clockwise  direction.  In  other  words,  the  orientation  of  the  trans¬ 
mitted  wave  will  decrease  if  the  thickness  of  the  chiral  slab  increases.  For  the 
case  of  a  wave  diffracted  by  a  chiral  grating,  the  same  argument  applies  for  each 
Floquet  mode.  When  the  wave  enters  gratings,  Bragg  diffraction  will  split  it  into 
Floquet  modes  propagating  in  several  directions,  the  orientation  of  each  mode 
will  rotate  in  the  clockwise  direction  for  the  case  of  right-hand  chirality.  However, 
it  should  be  noticed  that  for  some  particular  thickness,  the  orientation  of  some 
Floquet  mode  has  a  fast  transition  (Figs.  4-7(e)j  and  may  show  an  opposite  trend, 
such  as  in  Fig.  5(e)  where  the  orientation  of  Floquet  mode  0  increases  instead  of 
decreases  near  d  =  1.2A.  Nevertheless,  it  can  be  observed  that  whenever  there  is 
a  fast  transition  in  the  plots  of  the  orientation  angle  versus  depth,  the  correspond¬ 
ing  Floquet  mode  is  nearly  circularly  polarized,  for  example,  Floquet  mode  0  at  d 
near  1.2A  in  Figs.  4-7(f).  It  is  known  that  the  orientation  of  a  circularly  polarized 
wave  can  be  arbitrary  or  undefined.  Hence,  it  should  not  be  surprised  that  the 
orientation  exhibits  a  fast  transition  around  the  thickness  where  the  waves  are 
nearly  circularly  polarized.  Also  some  of  the  fast  jumps  shown  in  the  plots  of  ori¬ 
entation  angles  from  Odeg  to  180  deg  or  the  opposite  are  not  important  because 
Odeg  and  180  deg  orientations  represent  essentially  the  same  orientation. 

rv.  summary 

In  this  paper,  the  diffraction  of  waves  from  chiral  gratings  is  analyzed.  The 
problem  is  formulated  by  the  coupled  wave  method  which  gives  highly  convergent 
results  for  all  the  cases  studied  where  the  computed  total  diffraction  efficiency  is 
different  from  1  by  0.02%  for  the  worst  case.  It  is  shown  that  a  grating  of  chiral 
medium  is  capable  of  not  only  dividing  the  power  into  different  Floquet  modes, 
but  sometimes  also  making  both  dominant  modes  nearly  circularly  polarized  in 
the  cases  studied. 
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ABSTRACT 

We  present  an  inversion  algorithm  based  on  a  recently  developed  inversion  method  referred  to  as 
the  Renormalized  Source-Type  Integral  Equation  approach.  The  objective  of  this  method  is  to 
overcome  some  of  the  limitations  and  difficulties  of  the  iterative  Bom  technique.  It  recasts  the 
inversion,  which  is  nonlinear  in  nature,  in  terms  of  the  solution  of  a  set  of  linear  equations;  however, 
the  final  inversion  equation  is  still  nonlinear.  The  derived  inversion  equation  is  an  exact  equation  which 
sums  up  the  iterative  Neuman  (or  Bom)  series  in  a  closed  form  and;  thus,  is  a  valid  representation  even 
in  the  case  when  the  Bom  series  diverges;  hence,  the  name  Renormalized  Source-Type  Integral 
Equation  Approach. 


1.  INTRODUCTION 

Most  of  the  monochromatic  profile  inversion  methods  that  exist  in  literature  utilize  an  iterative  Bom 
approach  in  the  implementation  of  the  reconstruction  process^*®.  This  approach  does  not  work  well  at 
high  frequencies  or  for  the  transverse-magnetic  (TM)  polarization,  but  appears  to  be  useful  only  for  the 
low  frequency  transverse-electric  (TC)  case.  This  is  attributed  to  the  fact  that  at  high  frequencies  or  in 
the  case  of  the  TM  polarization,  the  inversion  problem  becomes  highly  nonlinear  rendering  perturbative 
approaches  inappropriate. 

In  this  paper  we  present  an  inversion  algorithm  based  on  a  recently  developed  inversion  method 
referred  to  as  the  Renormalized  Source-Type  Integral  Equation  approach^.lO.  The  objective  of  this 
method  is  to  overcome  some  of  the  limitations  and  difficulties  of  the  above-mentioned  iterative  Bom 
technique.  This  approach  can  be  summarized  in  three  simple  steps.  In  the  first  step,  the  induced 
currents  in  the  prob^  medium  are  determined  from  the  known  data  on  the  surface  of  medium  by 
use  of  the  source-type  integral  equation.  Next,  the  field  induced  in  the  medium  is  computed  from  the 
induced  currents.  Li  the  final  step,  the  medium  parameters  are  reconstructed  from  the  induced  currents 
and  fields  which  were  obtained  in  the  previous  two  steps.  This  method  recasts  the  inversion,  which  is 
nonlinear  in  nature,  in  terms  of  the  solution  of  a  set  of  linear  equations;  however,  the  final  inversion 
equation  is  still  nonlinear. 

There  are  a  number  of  apparent  advantages  to  this  scheme:  i)  it  is  applicable  to  TE,  TM  and  hybrid 
polarizations,  ii)  there  is  no  explicit  frequency  limitations,  iii)  by  pre-storing  the  elements  of  the 
inversion  operator  (which  depend  only  on  the  background  me^um  and  are  independent  of  the  unknown 
profile),  the  method  does  not  require  the  solution  to  the  full  forward  problem  repeatedly  as  in  the  case 
of  the  Distorted  Bom  approach  and  is  therefore  faster  to  implement,  iv)  it  can  be  generalized  to  two- 
dimensional  cases  where  the  previous  property  becomes  an  attractive  feature  of  the  inversion,  v)  it 
allows,  in  some  cases,  the  rigorous  study  of  the  degree  of  nonuniqueness  involved  in  the  inversion,  a 
point  of  great  importance  in  tiie  design  problem,  and  vi)  it  provides  an  inversion  equation  which  has  an 
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explicit  dependence  on  the  unknowns  to  be  obtained  in  the  inversion.  This  explicit  dependence  allows 
one  to  compute  the  derivatives  of  the  response  with  respect  to  the  unknowns  in  a  closed  form.. 
Associated  with  this  approach,  however,  is  the  need  to  incorporate  a  priori  information,  such  as  the 
type  of  likely  profiles,  that  are  derived  from  means  other  than  those  supplied  by  measurement.  The 
judicious  use  of  such  information  allows  us  to  overcome  the  non-uniqueness  which  is  characteristic  of 
all  inversion  problems. 

Another  advantage  of  this  approach  is  that  the  inversion  equation  obtained  is  an  exact  equation.  If 
approximations  are  to  be  made  they  will  occur  in  the  numerical  methods  employed  in  the  solution  of  this 
inversion  equation.  Thus,  the  major  virtue  of  our  method  is  that  the  approximations  to  be  used  can  be 
under  good  control  and  may  be  improved  upon  to  an  arbitrary  degree.  This  is  to  be  distinguished  from 
other  methods,  e.g..  Bom  inversion,  ray  tracing,  etc.,  which  are  inherently  approximate,  often  unstable 
and  hard  to  improve  when  they  give  inaccurate  results. 

In  this  paper,  the  inversion  scheme  is  demonstrated  for  the  cases  where  the  exciting  source  is  of  the 
dipole-type  that  generates  only  TE  or  TM  polarized  waves.  The  data  required  for  inversion  are  the 
tangential  components  of  the  electric  or  magnetic  fields  at  various  locations  on  the  probed  medium's 
surface  performed  at  a  number  of  frequencies. 

2.  FORMULATION 

We  limit  our  analysis  to  the  problem  of  profile  inversion  of  a  one-dimensional  layered  medium 
whose  permittivity  can  vap'  as  a  function  of  depth  in  a  planar  stratified  medium  or  as  a  function  of 
radial  distance  in  a  cylindrically  layered  medium.  For  the  sake  of  simplifying  the  analysis,  we  restrict 
ourselves  to  the  planar  stratified  case  where  the  medium  is  assumed  to  be  z-stratified,  i.e.,  the 
permittivity  e  and  the  conductivity  o  of  the  medium  are  functions  of  z  only.  The  geometrical 
configuration  of  the  problem  is  depicted  in  figure  1. 
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Figure  1.  Geometrical  configuration  of  the  inversion  problem. 

The  source  is  assumed  to  be  localized,  and  without  any  loss  of  generality,  the  source  characteristics 
can  be  chosen  such  that  they  excite  either  a  purely  transverse  electric  (TE)  or  transverse  magnetic  (TM) 
type  of  field.  For  the  sake  of  further  simplifying  the  analysis,  we  assume  that  the  exciting  TE  source  is 
an  electric  line  current,  situated  at  x  =  0,  z  =  Zs  and  extending  from  y  =  -  oo  to  y  =  -f  oo,  in  this  case, 
the  nonvanishing  field  components  are  Hx,  Hz  and  Ey  which  are  functions  of  x  and  z  only.  For  the  TM 
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case,  we  assume  the  source  to  be  a  magnetic  line  current  and,  in  this  case,  the  nonvanishing  field 
components  are  Ex,  Ez  and  Hy. 

In  formulating  the  problem  using  the  source-type  integral  equation,  one  assumes  that  the  unknown 
inhomogeneous  slab  (extending  from  z  =  0  to  z  =  L)  is  embedd^  inside  a  known  background  medium 
(ejj(z))  which  is  not  necessarily  homogeneous.  The  scattering  from  this  inhomogeneous  slab  is  then 
assumed  to  be  produced  by  a  volume  distribution  of  currents  induced  inside  the  slab.  The  field  can  thus 
be  decompose  into  a  superposition  of  two  parts.  The  first  part  (Eyjj  for  TE  and  Hy^,  for  TM)  is 
produced  by  the  exciting  source  in  the  background  medium.  The  second  (Ey^  for  TE  and  Hyj  for  TM) 
represents  the  field  scattered  by  the  inhomogeneous  slab  and  is  given  by  the  field  produced  by  the 
volume  currents  induced  inside  the  slab.  Thus,  for  the  TE  polarization  we  have: 

Ey  =  Eyb  Ey,  (1) 

V^Eyb-»-k^e^,Eyb  =  0  (2a) 

V^Ey.  -h  ko^  Eb  Ey,  =  -ko^  5e  Ey  (2b) 


and  for  the  TM  polarization,  we  have: 


Hy  =  Hyb  +  Hy, 


.  £yf  dHyb  ^ 


V^H,.  -  +  ko^  e,  =  -ko^  5e  H,  + 


e;_ 

e 


iL 

Eb 


aHy 

dz 


(3) 

(4a) 

(4b) 


where 


5e  =  e(z)  -  eb(z)  (5) 

Using  the  Green's  function  formulation,  eqs.  (2b)  and  (4b)  can  be  cast  in  the  following  integral  form: 


Ey,(x,  z)  =  k^  J  dz'  5e  J  dx'  G  (x,  z,  x',  z')  Ey(x',  z') 
0  -- 


(6) 


for  the  TE  polarization.  For  the  TM  polarization,  we  get 


Hy,(x,z)  =  J  dz'  5e  J  dx'  G^(x,  z,  x',  z')  Hy(x',z') 


-  J  dz' 


ll  _iL 

e  Eb 


J  dx'  G^(x,  z,  x',  z')  —  Hy(x',  z') 


(7) 
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where  Gjg  and  Gjj^  are  the  Green's  functions  for  the  TE  and  TM  polarizations,  respectively,  and  arc 
given  as  the  solution  of: 


Cb  =  -5(x-x0  5(z-z')  (8a) 

The  source-type  integral  equation  is  then  solved  by  applying  successive  iterations  where  a 
Distorted-Bom  approximation  is  applied  at  each  iteration  step.  For  the  iteration  scheme  to  converge, 
lEyjl  and  IHy^l  have  to  be  small.  It  is  clear  that  the  term  on  the  right-hand  side  of  eq.  (6)  and  the  first 
term  on  the  right-hand  side  of  eq.  (7)  can  be  made  arbitrarily  small  by  either  lowering  the  frequency  of 
operation  Icq  or  making  I5el  small.  However,  the  second  term  on  the  right-hand  side  of  eq.  (7)  can  only 
be  made  small  by  making  I5el  small,  i.e.,  if  t^(z)  is  sufficiently  close  to  £(z).  Thus,  one  can  conclude 
that  by  choosing  a  frequency  of  operation  that  is  sufficiently  low,  the  iteration  can  be  made  to  converge 
in  the  case  of  the  TE  polarization  even  for  fairly  large  values  of  I5eI  ,  whereas  in  the  case  of  the  TM 
polarization,  the  only  way  the  iteration  will  converge  is  when  the  initial  guessed  profile  e^fz)  is  chosen 

to  be  very  close  to  the  true  profile  e(z).  It  is  interesting  to  note  that  the  presence  of  the  second  term  on 
the  right-hand  side  of  eq.  (7),  causes  the  TM  measurement  to  be  iTX>re  sensitive  to  abrupt  changes  in  the 
permittivity  profile.  Thus,  the  TM  measurement  is  a  better  candidate  than  the  TE  case  for  imaging 
applications  where  one  can  delineate  boundaries. 

Another  difference  between  the  TE  and  TM  polarizations  lies  in  the  degree  of  nonlinearity  associated 
with  each  case.  It  is  clear  that  eqs.  (6)  and  (7)  are  nonlinear  integral  equations  in  the  unknown  profile 

e(z)  due  to  the  nonlinear  dependence  of  the  field  (Ey  or  Hy)  on  e(z).  Using  a  Bom-type  approximation 
in  carrying  out  the  iteration  is  in  effect  linearizing  tne  problem.  This  linearization  is  expected  to  work 
better  in  the  case  of  the  TE  polarization  than  in  ^e  TM  case  because  the  TM  case  is  highly  nonlinear 
compared  to  the  TE  case.  TTiis  can  be  seen  by  transforming  eqs.  (2)  and  (4)  into  a  canonical  form  for 
comparison.  This  is  done  by  first  carrying  out  a  Fourier  transform  along  the  x-direction  (assuming  k^ 
to  be  the  transform  variable)  and  then  introducing  the  field  transformations 


u(k„  z)  =  E(kj,  z)  , 
u(k,, z)  =  H(k,, z) , 

This  yields  the  Schroedinger-typ)e  equation  in  the  field  variable  u 

9^  2 

— -  +  Y  u  =  -qu 
az2 

where  y  is  given  by 

Y^  =  k^  Eb  -  ki 


for  TE 

(9a) 

for  TM 

(9b) 

(10) 

(11) 

SPI£  Vo/  1351  Digital /mage  Synthesis  and  Inverse  Optics  (1 990)  /  101 


and 


q  =  k^  (E-Eb) , 


q  =  k^  (E-Eb)-  — 
4 


e  ^ 
E' 


-T-- 

2  E 


for  TE 

(12a) 

for  TM 

(12b) 

Note  that  the  right-hand  side  of  cq.  (10)  can  be  interpreted  as  a  source  term  resulting  from  the 
medium  inhomogeneities  in  the  region  z  >  0.  From  eq.  (12),  it  is  thus  clear  that  for  the  TE  polarization, 
e(z)  is  direcdy  calculable  from  q;  however,  it  is  necessary  to  solve  a  second-order  nonlinear  differential 
equation  to  obtain  e(z)  for  the  case.  TTius,  the  TM  case  is  relatively  nonlinear  compared  to  the  TE 
case. 


3.  THE  RENORMALT7.FD  SOT  TRCE-TYPE  INTEGRAL  EOT  JATION  APPROACH 

The  Renormalized  Source-Type  Integral  Equation  approach^**®  can  be  summarized  in  the  following 
three  simple  steps: 

Step  1.  The  Inverse  Source  Problem:  From  the  collected  data  on  the  surface  of  the  probed  medium 
(r  E  Sm),  the  currents  induced  inside  the  medium,  J(r),  are  solved  for  from  the  following  equation: 

/drGb(T;-r)T(r)  =  E(T)-Eb(T)|^^s_  (i3) 

Step  2.  The  Forward  Problem:  From  the  knowledge  of  the  currents  induced  inside  the  slab  (obtained 
from  step  1),  the  electric  field,  E(r),  can  be  calculated  everywhere  (inside  as  well  as  outside  the  slab) 
from  the  following  equation; 

E(T)  =  Eb(T)  +  \dr  Gb(  r,r)-T(r)  (14) 

Step  3.  The  Inversion  Equation:  From  the  knowledge  of  the  currents  induced  inside  the  slab 
(obtained  from  step  1)  and  the  electric  field  inside  the  slab  (computed  from  step  2),  the  unknown 
difference  profile  Q(r)  =  (o^pQ[E(r)  -  eJ  is  obtained  from  the  following  constitutive  relation: 

Q(r)E(T)=T(T)  (15) 


Thus,  this  method  recasts  the  inversion,  which  is  nonlinear  in  nature,  in  terms  of  the  solution  of  a 
set  of  linear  equations;  however,  the  final  inversion  equation  is  still  nonlinear.  In  this  paper,  the 
inversion  scheme  is  demonstrated  for  the  cases  where  the  exciting  source  is  of  the  dipole  type  that 
generates  only  TE  or  TM  polarized  waves.  The  data  required  for  inversion  are  the  tangential 
components  of  the  electric  and/or  magnetic  fields  at  various  locations  on  the  surface  of  the  probed 
medium  performed  at  a  number  of  frequencies.  In  the  next  section  this  inversion  approach  is  tested  on 
a  one-dimensional  profile  where  data  are  generated  by  exciting  the  medium  once  by  a  horizontal 
magnetic  dipole  (exciting  TE-polarized  waves)  and  a  second  time  by  a  horizontal  electric  dipole 
(exciting  TM-polarized  waves). 
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In  this  section  we  apply  the  inversion  algorithm  described  in  the  previous  section  to  the  problem 
depicted  in  figure  2. 
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Figure  2.  Permittivity  profile  used  in  the  inversion 


The  data  used  in  the  inversion  are  generated  by  exciting  the  medium  once  by  a  horizontal  magnetic 
dipole  (exciting  TE-polarized  waves)  and  a  second  time  by  a  horizontal  electric  dipole  (exciting  TM- 
polarized  waves).  The  data  are  collected  at  ten  receiver  locations  with  the  appropriate  type  and 
orientation.  A  single  electric  field  component  (E,j,)  is  measured  for  the  TE  case  whereas  two  electric 
field  components  (Ez  and  Ep)  are  measured  for  the  TM  case.  Note  that  the  dimensions  of  the  slab  is 
normalized  to  the  space  wavelength. 

Figure  3a  shows  the  result  of  inversion  using  only  TE-polarized  data  whereas  figure  3b  shows  the 
result  of  inversion  using  only  TM-polarized  data.  From  figure  3a,  it  is  clear  that  the  reconstructed 
profile  employing  TE  data  is  a  filtered  (or  a  smoothed)  version  of  the  true  profile  and  hence  TE- 
polarized  data  offer  poorer  resolution  in  delineating  abrupt  changes  in  the  permittivity  profile.  On  the 
other  hand  TM-polarized  data  offer  a  higher  resolution,  however,  the  invert^  profile  exhibits  nonlinear 
effects  depicted  by  the  presence  of  hom-type  features  in  the  reconstructed  profile.  From  figure  3,  it  is 
also  clear  that  TE-inversion  offers  higher  accuracy  than  TM  in  retrieving  the  permittivity  values  around 
the  peak  of  the  profile.  It  is  also  interesting  to  note  that  the  TM-inversion  has  a  shallow  depth  of 
investigation  which  is  manifested  by  the  lower  resolution  at  the  further  edge  of  the  profile  compared  to 
the  higher  resolution  at  the  near  edge.  Finally,  the  TE-based  inversion  converged  after  8  iterations 
whereas  the  TM-based  inversion  converged  after  33  iterations.  Hence,  faster  convergence  can  be 
achieved  by  employing  TE-polarized  data. 
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Figure  3b.  Inversion  result  employing  TM-polarized  CW  data. 
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Figure  4  shows  the  result  of  inversion  combining  both  TE  and  TM  polarized  data  where  improved 
resolution  and  accuracy  has  been  achieved. 


Figure  4.  Inversion  result  employing  the  combined  TE  and  TM  polarized  CW  data. 

5.  CONCLUSIONS 

In  this  paper  we  have  presented  a  new  and  exact  inversion  equation  which  is  applicable  to  data  of 
TE,  TM  or  hybrid  polarization.  If  approximations  are  to  be  made  they  will  occur  in  the  numerical 
methods  employed  in  the  solution  of  this  inversion  equation.  This  is  to  be  distinguished  from  other 
methods,  e.g..  Bom  inversion,  ray  tracing,  etc.,  which  are  inherently  approximate  and  are  hard  to 
improve  when  they  give  inaccurate  results.  The  inversion  equation  is  a  closed  form  summation  of  the 
iterative  (or  perturbative)  Neumann  (or  Bom)  series  of  the  Source-Type  Integral  Equation.  This 
equation  is  a  valid  representation  even  in  the  case  when  the  Bom  series  diverges. 
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Because  the  effects  of  diffraction  during  proximity-print  x-ray 
lithography  are  of  critical  importance,  a  number  of  previous 
researchers  have  attempted  to  calculate  the  diffraction  patterns  and 
minimum  achievable  feature  sizes  as  a  function  of  wavelength  and 
gap.  Work  to  date  has  assumed  that  scalar  diffraction  theory  is 
applicable — as  calculated,  for  example,  by  the  Rayleigh-Sommerfeld 
formulation — and  that  Kirchhoff  boundary  conditions  can  be  applied. 
Kirchhoff  boundary  conditions  assume  that  the  fields  (amplitude  and 
phase)  are  constant  in  the  open  regions  between  absorbers,  and  a 
different  constant  in  regions  just  under  the  absorbers  (i.e.,  that  there  are 
no  fringing  fields).  An  x-ray  absorber  is,  however,  best  described  as  a 
lossy  dielectric  that  is  tens  or  hundreds  of  wavelengths  tall,  and  hence 
Kirchhoff  boundary  conditions  are  unsuitable.  In  this  report  we  use 
two  numerical  techniques  to  calculate  (on  a  Cray  2  supercomputer) 
accurate  diffracted  fields  from  gold  absorbers  for  two  cases:  a  30  nm- 
wide  line  at  X  =  4.5  nm,  and  a  100  nm-wide  line  at  X  =  1.3  nm.  We 
show  that  the  use  of  Kirchhoff  boundary  conditions  introduces 
unphysically  high  spatial  frequencies  into  the  diffracted  fields.  The 
suppression  of  these  frequencies — which  occurs  naturally  without  the 
need  to  introduce  an  extended  source  or  broad  spectrum — improves 
exposure  latitude  for  mask  features  near  0.1  ^im  and  below. 
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I.  INTRODUCTION 


The  limits  and  practicality  of  proximity-print  x-ray  lithography  continues  to 
be  a  topic  of  discussion  and  debate.  Of  particular  concern  are  the  limits  of 
resolution  imposed  by  the  effects  of  diffraction.  Because  the  mechanical 
limits  imposed  on  the  mask-substrate  gap  during  the  volume  manufacturing 
of  ULSI  circuits  are  not  presently  certain,  it  remains  prudent  to  ask:  what  is 
the  minimum  practical  feature  size  that  can  be  printed  at  a  given  gap? 

In  order  to  resolve  these  issues,  a  series  of  papers  have  appeared  in  the 
literature  which  present  theoretical  calculations  for  the  diffraction  of  x  rays 
from  mask  absorbers.  Early  papers  simply  considered  absorption  [1-3],  but 
later  papers  also  included  the  effects  of  phase-shift  [4-8].  Most  recently, 
authors  have  included  the  effects  of  source  spatial  and  temporal  coherence, 
and  have  generated  exposure  "trees"  which  allow  the  determination  of 
exposure  latitude  versus  gap  for  various  types  of  mask  features  [9-15]. 

The  method  most  commonly  used  to  calculate  the  diffraction  pattern  is  to 
apply  a  Fresnel-Kirchhoff  or  Rayleigh-Sommerfeld  diffraction  integral  [16,17], 
or  a  more  sophisticated  formulation  (Hopkin's  formula)  which  takes  into 
accovmt  source  partial  coherence  [18].  These  calculations  can  be  carried  out 
either  in  the  spatial  or  the  spatial-frequency  domain.  In  any  of  these  cases, 
approximate  boundary  conditions  known  as  Kirchhoff  boundary  conditions 
(KBC)  are  generally  applied.  KBC  assume  that  the  field  (amplitude  and  phase) 
is  constant  in  the  open  regions  between  absorbers,  and  also  constant  (but 
attenuated  and  phase  shifted)  in  regions  just  under  the  absorbers — in  other 
words,  that  there  are  no  fringing  fields  at  the  boundary  between  the  two 
regions. 

In  general,  KBC  apply  when  the  wavelength,  X,  is  much  smaller  than  the 
lateral  size,  d,  of  the  feature  being  printed — which  is  the  case  in  most  x-ray 
lithography  (e.g.,  d  =  0.1  pm,  X  =  1  nm).  However,  what  is  not  generally 
recognized  is  that  KBC  will  not  necessarily  apply  when  the  absorbers  are  lossy 
dielectrics  that  are  tens  or  hundreds  of  wavelengths  tall. 

Work  to  date  concerning  the  printability  of  0.25  pm  features  is  probably 
reasonably  accurate.  However,  in  this  report  we  show  that  the  assumption  of 
KBC  for  features  near  0.1  pm  and  below  is  not  tenable.  In  particular,  we  show 
that  the  suppression  of  the  undesirable  high-spatial-frequency  components — 
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which  some  authors  note  an  extended  source  and/or  a  broad  spectrum 
achieves — occurs  natvu-ally  in  the  absorber  due  to  the  "lossy  dielectric"  effect. 

IL  CALCULATIONS 

We  used  two  different  methods  to  calculate  the  diffraction  from  gold 
absorbers:  the  Method  of  Moments  (MoM)  and  the  Finite  Difference-Time 
Domain  (FD-TD)  method.  (We  used  two  algorithms  in  order  to  check  for 
consistency.)  The  only  approximations  inherent  in  these  methods  are  in  the 
discretization  of  the  object  space  and  Maxwell's  equations.  Because  the 
discretization  can  performed  on  a  scale  that  is  small  compared  to  the 
wavelength,  and  furthermore  the  discretization  scale  can  be  reduced  until 
convergence  is  achieved,  accuracy  is  assured. 

Of  these  two  techniques,  the  MoM  is  typically  faster  and  uses  less  memory 
for  single-frequency  calculations.  On  the  other  hand,  the  FD-TD  method  is 
simpler  to  code  and  therefore  less  likely  to  have  errors.  In  practice  we  ran 
small  test  cases  using  both  methods  and  then  increased  the  spatial  resolution 
(reduced  discretization  scale)  until  both  methods  converged  to  the  same 
solution.  Then  the  computationally-intensive  cases  reported  here  were 
calculated  with  the  MoM. 

A.  METHOD  OF  MOMENTS 

The  Method  of  Moments  (MoM)  is  a  numerical  technique  useful  in  the 
solution  of  steady-state  electromagnetic  wave  scattering  and  radiation 
problems  [19,20].  The  method  calculates  steady-state  fields  on  the  surface  of  a 
closed  dielectric  object  in  free  space,  which  is  impinged  upon  by  a  known 
exciting  wave.  The  surface  of  the  object  is  broken  up  into  small  patches 
which  are  small  compared  to  the  wavelength.  The  surface  currents  at  each 
patch  and  thiw  the  tangential  surface  fields  are  then  calculated.  Computation 
time  was  up  to  two  hours  on  a  Cray  2.  Once  the  fields  are  knovm  on  the 
boimdary  of  the  object  they  can  readily  be  propagated  to  any  desired  point  or 
plane  [16,17]. 
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B.  FINITE  DIFFERENCE-TIME  DOMAIN  METHOD 


The  Finite  Difference-Time  Domain  (FD-TD)  method  is  a  numerical 
technique  useful  in  the  solution  of  time-dependent  electromagnetic  wave 
scattering  and  radiation  problems  [21-24].  The  method  involves  the 
formation  of  a  computational  domain  which  encompasses  the  object  of 
interest  and  is  typically  several  times  larger.  The  entire  domain  inside  the 
boundary — including  the  object — is  discretized  on  a  rectangular  grid.  The 
spacing  between  adjacent  nodes  on  the  grid  is  small  compared  to  the 
wavelength.  A  clock  is  started  and  incremented  in  time  steps  that  are  small 
compared  to  the  light-travel  time  between  adjacent  nodes.  Then  discretized 
forms  of  Maxwell's  equations  are  used  to  calculate  the  fields  at  each  node 
from  the  fields  at  nearby  nodes  which  were  in  effect  at  the  previous  time  step. 
Absorbing  boundary  conditions  are  imposed  at  the  edges  of  the 
computational  domain  in  order  to  simulate  unbounded  space.  Also,  a 
boundary  condition  is  typically  imposed  on  a  surface  surrounding  the  object 
to  simulate  an  incoming  plane  wave.  The  result  can  be  displayed  as  a  video 
image  of  the  fields  inside  the  domain.  Computation  time  was  up  to  one  hour 
on  a  Cray  2.  In  practice  the  calculation  is  run  until  steady  state  is  achieved, 
and  the  fields  at  nodes  along  a  line  just  under  the  absorber  are  saved  for 
comparison  with  the  MoM  results. 

III.  RESULTS 

Calculations  were  performed  using  MoM  and  FD-TD  on  single  gold 
parallelepiped  absorbers,  infinite  in  length  and  rectangular  in  cross-section, 
which  were  impinged  upon  by  a  monochromatic  plane  wave.  The  electric- 
field  polarization  was  used  (E-field  perpendicular  to  the  page).  (The 
magnetic-field  polarization  yielded  results  similar  to  the  electric-field 
polarization.)  These  were  compared  to  the  results  of  a  Rayleigh-Spmmerfeld- 
Kirchhoff  (RSK)  calculation  [16].  We  considered  two  cases  (see  Table  I),  which 
were  selected  in  order  to  explore  a  range  of  spatial  frequencies:  Case  1,  which 
is  a  30x100  nm  absorber  with  the  4.5  nm  (Ck)  x  ray,  and  Case  2,  which  is  a 
100x250  nm  absorber  with  the  1.3  nm  (Cul)  x  ray.  Note  that  the  attenuation  of 
the  absorber  in  both  cases  is  roughly  equivalent.  (~12  dB). 
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A.  CASE  1  a  =  4.5  lun) 


Figs.  1  (a)  to  (d)  show  the  resulting  intensities  or  irradiance  distributions 
for  both  the  MoM  and  RSK  methods  calculated  for  a  from  0-1.5  (gap  G  from  0- 
0.3  pm).  Here  a  is  a  dimensionless  gap  given  by  a  =  GX/W^,  where  G  is  the 
gap,  X  is  the  wavelength,  and  W  is  the  linewidth.  Note  the  suppression  of 
high  spatial  frequencies  in  the  MoM  calculation,  and  the  "fuzzy  edge"  of  a  few 
tens  of  nanometers  in  extent.  Figs.  2  (a)  and  (b)  show  the  beneficial  effects  of 
the  suppression  of  high  spatial  frequencies  on  exposure  latitude  in  the  form 
of  exposure  "trees"  [9-12]  which  plot  ±10  %  linewidth  contours  versus  a.  Here 
we  have  used  a  line  bias  of  33  %  (40  nm  resist  line).  (A  line  bias  is  the  use  of  a 
smaller-than-desired  feature  size  in  the  mask  than  on  the  wafer  in  order  to 
compensate  for  diffractive  broadening.)  An  enlightening  way  to  view  the 
suppression  of  high  spatial  frequencies  is  to  plot  the  intensity  of  the  waves  at 
a  =  0  (gap  G  =  0)  as  a  function  of  spatial  frequency,  as  shown  in  Fig.  3.  Note 
that  at  low  spatial  frequencies  the  RSK  and  MoM  calculations  agree,  but  the 
MoM  calculation  "rolls  off"  at  around  0.05  nm*'. 

B.  CASE  2  (X  =  1.3nm) 

The  results  for  the  intensities  versus  a  for  Case  2  are  shown  in  Figs.  4  (a)  to 
(d).  Note  that  even  though  the  wavelength  is  much  smaller,  the  "fuzzy  edge" 
length  is  roughly  the  same — a  few  tens  of  nanometers.  This  may  be  due  to 
the  smaller  wavelength  (3.3x  smaller)  being  partially  compensated  by  a  taller 
absorber  (2.5x  taller).  The  exposure  trees  are  shown  in  Figs.  5  (a)  and  (b).  Here 
we  have  used  a  line  bias  of  50%  (150  nm  resist  line).  The  intensity  versus 
spatial  frequency  is  shown  in  Fig.  6.  Note  that  the  roll-off  in  this  case  is  still 
around  0.05  nm*^,  but  that  this  represents  a  higher  spatial  frequency  relative 
to  the  information  content  in  the  larger-line/smaller-wavelength  case. 

IV.  CONCLUSION 

We  have  shown  that  the  use  of  Kirchhoff  boundary  conditions  introduces 
unphysically  high  spatial  frequencies  into  the  diffracted  fields.  The  natural 
suppression  of  these  frequencies  by  the  electromagnetic  properties  of  x-ray 
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absorbers  tremendously  improves  exposure  latitude  for  mask  features  near 
0.1  pm  and  below. 
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Figure  Captions 

FIG.  1.  Intensity  vs.  lateral  position  for  Case  1.  (a)  a  =  0  (gap  =  0),  (b)  a  =  0.5 
(gap  =  0.1  pm),  (c)  a  =  1.0  (gap  =  0.2  pm),  (d)  a  =  1.5  (gap  =  0.3  pm). 

FIG.  2-  Exposure  trees  vs.  a  (dimensionless  gap)  for  Case  1.  (Gaps  range  from 
0-0.3  pm.)  The  line  is  biased  33%  (40  nm  resist  line),  (a)  MoM 
solution,  (b)  Rayleigh-Sommerfeld-Kirchhoff  approximation. 

FIG.  3  Intensity  vs.  spatial  frequency  for  Case  1.  Note  roll-off  at  0.05  nm*'. 

FIG.  4  Intensity  vs.  lateral  position  for  Case  2.  (a)  a  =  0  (gap  =  0),  (b)  a  =  0.5 
(gap  =  3.7  pm),  (c)  a  =  1.0  (gap  =  7.5  pm),  (d)  a  =  1.5  (gap  =  11.2  pm). 

FIG.  5  Exposure  trees  vs.  a  (dimensionless  gap)  for  Case  2.  (Gaps  range  from 
0-11.2  pm.)  The  line  is  biased  50%  (150  nm  resist  line),  (a)  MoM 
solution,  (b)  Rayleigh-Sommerfeld-Kirchhoff  approximation. 

FIG.  6  Intensity  vs.  spatial  frequency  for  Case  2.  Note  roll-off  at  0.05  nm*'. 
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TABLE  I.  Computational  cases  considered. 


Case  1  Case  2 

Wavelength  X  (run)  4.48  1.334 

Refractive  Index  [25]  1  -  7.54x10-3  +  jl.04xl0-2  1  -  2.31x10-3  +  jl.:9xl0-3 

Width  W  (nm)  30  100 

Height  (nm)  100  250 

Transmission  0.0541  0.0607 

Phase  Shift  (rad)  1.058  2.72 

Patch  Size  (run)  -0.64  -0.19 
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